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KXKCUTIVK  SUMMARY 


This  document  presents  an  annual  report  to  the  Office  of  Naval  Research  for  a  research 
program  entitled  “A  Theoretical  Search  for  Supervelocity  Semiconductors.”  This  program  has 
been  funded  by  ONR  since  1974  in  the  Department  of  Electrical  and  Computer  Engineering  at 
N.C  State  University.  The  research  has  resulted  in  more  than  75  refereed  publications  and 
numerous  conference  presentations  from  its  inception.  Major  contributions  to  the  field  of  hot 
electron  iranspoit  and  semiconductor  device  modeling  have  been  achieved,  new  computational 
methods  have  been  developed  (e.g.,  path  integral  Monte  Carlo  techniques),  and  the  wc.rk  has 
helped  stimulate  commercial -ventures  in  the  applications  of  quaternary  semiconductor materials  to 
electronic  and  optical  devices.  In  addition,  there  have  been  twenty-six  Ph.D.  and  M.S.  students 
who  have  received  degrees  at  N.C.  State  University  with  research  support  from  inis  contract. 
Three  visiting  faculty  members  from  Japan  came  to  the  University  to  work  with  the  faculty 
investigators  supported  under  this  ONR  contract  during  the  1979-1983  time  period.  A  visiting 
professor  from  die  French  CNRS  Micros  true  lures  and  Microelectronics  Laboratory  in  Bagneux 
(near  Paris)  spent  a  sabbatical  year  at  N.C.  State  during  1988-89.  and  he  devoted  full-time 
working  on  this  program  at  no  cost  to  ONR.  During  the  current  funding  period,  a  visiting  scholar 
from  China  is  a  member  of  our  research  group  working  on  projects  which  directly  impact  this 
ONR  program.  Furthermore,  we  expect  to  add  a  new  Visiting  Research  Assistant  Professor  this 
year  to  enhance  our  efforts  in  band  structure  calculations  for  quantum  transport  and  mesoscopic 
devices. 


The  initial  phase >  of  this  work  centered  around  development  of  Monte  Carlo  simulation 
techniques  which  allow  the  study  of  detailed  physics  of  hot  electron  transport  :n  a  variety  of 
compound  semiconductor  materials.  The  original  emphases  were  condensed  .vuh  electron!, 
materials  phenomena  Later  work  considered  the  ii.ili/ation  of  these  materials  i»  realistic  ilcv  ».v 
structures  where  physical  boundary  conditions  mu  a  be  imposed  on  the  vutnei  transport  Mote 
recently,  the  work  ha*  locw.ednn  the  domaui  ofulira  sir.il!  materi.:i.%  and  device  phenomena  v.h  ”. 
microscopic  non  local  tianstcnl  CmCIsmuIi  ,i-  veloeitv  overshoot,  ballisi:.  a:n.  nearly  balhstu 


transport,  and  quantum  transport  become  important  or  dominant.  During  the  past  four  years  we 
have  published  new  research  results  on  applications  of  the  Feynman  “integral  over  paths”  approach 
to  quantum  transport  as  well  as  the  study  of  hot  electron  effects  in  new  device  structures,  such  as 
the  hot  electron  spectrometer,  heterojunction  bipolar  transistor,  small  dimension  metal- 
semiconductor-metal  photodetectors,  and  delta-doped  high  electron  mobility  transistors.  Also,  we 
have  been  exploring  some  new  approaches  to  device  modeling  which  combine  the  Monte  Carlo 
method  with  the  method  of  moments  of  the  Boltzmann  transport  equation  (hydrodynamic  transport 
model)  for  studying  specific  device  structures,  such  as  small-dimension  n+-n-n+  majority  carrier 
devices  and  the  high  electron  mobility  transistor.  During  the  past  year  we  have  incorporated 
quantum  correction  terms  into  the  hydrodynamic  model  and  applied  this  model  to  resonant 
tunneling  structures.  Quite  recently,  we  have  applied  a  new  ansatz  distribution  function  as  a 
constitutive  relation  to  close  the  moment  equations  in  the  hydrodynamic  transport  model.  Initial 
results  of  this  approach  have  been  physically  satisfying  and  computationally  promising.  A 
comprehensive  formulation  of  this  new  model  will  be  published  in  the  Physical  Review  in  January, 
1992.  New  work  on  the  theory  of  optical  phonon  modes  in  heterostructures  has  demonstrated  the 
importance  of  the  effects  of  reduced  dimensionality  on  transport  dominated  by  the  LO  phonon 
interaction  near  heterointerfaces.  As  a  result  of  this  work,  a  new  concept  for  velocity  enhancement 
in  pseudomorphic  device  structures  is  proposed  for  the  new  research  program.  Finally,  we  have 
developed  a  concept  for  a  new  device  structure,  which  we  call  the  real  space  transfer  logic 
transistor  (RSTLT).  Preliminary  results  for  this  dev,r  e  predict  a  switching  speed  less  than  3ps. 
These  results  will  be  presented  at  the  International  Electron  Devices  Meeting  (IEDM)  in  December. 
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1.0  Introduction 


In  October  1974,  the  Office  of  Naval  Research  initiated  sponsorship  of  a  basic  research 
program  in  the  Department  of  Electrical  and  Computer  Engineering  (ECE)  at  North  Carolina  State 
University.  The  general  goal  of  this  research  program  has  been  and  continues  to  be  the 
investigation  of  high-speed  carrier  transport  in  III-V  compound  semiconductors,  III-V  alloy 
materials,  other  advanced  electronic  and  optical  materials,  and  novel  device  structures  which  utilize 
these  materials.  Four  faculty  members  at  N.C.  State  have  been  primarily  involved  in  and 
supported  by  this  project.  They  are  M.A.  Littlejohn,  J.R.  Hauser  and  T.H.  Glisson,  presently 
Professors  in  the  ECE  Department,  and  Dr.  K.W.-Kim  who  joined  the  ECE  faculty  in  August 
1988.  Dr.  Kim  has  served  as  a  PI  on  the  effort  since  January  1989  and  is  now  a  major  contributor 
to  these  efforts.  This  research  program  has  made  significant  contributions  to  the  understanding 
and  knowledge  base  of  hot  electron  transport  in  materials  and  devices.  It  has  provided  scientific 
guidance  to  the  U  S.  Navy  in  the  formulation  of  a  part  of  its  basic  and  applied  research  program. 
Numerical  concepts  developed  under  this  project  have  been  transferred  to  other  Universities, 
including  the  University  of  Illinois,  and  we  continue  these  intcrinstitutional  collaborations.  During 
the  past  year,  faculty  members  and  graduate  students  have  spent  time  at  the  National  Center  for 
Computational  Electronics  at  the  University  of  Illinois  (with  Professor  K.  Hess).  Also, 
collaborations  with  Dr.  ILL.  Grubin  of  Scientific  Research  Associates,  Inc.,  Dr.  M.A.  Stroscio 
and  Dr.  G.J.  lafrate  of  the  U.S.  Army  Research  Office  in  the  area  of  quantum  transport  in 
semiconductor  devices  have  been  strengthened.  Dr.  lafrate,  Dr.  Stroscio,  and  Dr.  Grubin  are 
Adjunct  Professors  in  the  Electrical  and  Computer  Engineering  department  at  North  Carolina  State 
University.  In  addition  to  those  directions,  our  research  results  have  helped  stimulate  commercial 
ventures,  particularly  in  the  development  of  GalnAsP-based  materials  and  devices.  To-date,  this 
program  has  resulted  in  68  refereed  publications  in  the  literature,  6  additional  manuscripts  are 
currently  in  press,  and  5  manuscripts  have  been  submitted  or  are  in  preparation  for  submission  to 
the  technical  literature.  A  listing  of  these  publications  is  given  in  Appendix  A.  In  addition, 
numerous  invited  talks  and  presentations  have  been  given  at  conferences  and  workshops 
throughout  the  United  States  and  in  other  countries.  The  program  has  contributed  significantly  to 
the  educational  program  at  N.C.  State  University  with  more  than  thirty  Ph.D.,  M.S.,  and 
undergraduate  students  having  received  support  unde,  this  ONR  contract.  Currently,  three  Ph.D. 
students  are  working  toward  their  degrees  on  this  project.  Two  of  these  Ph.D.  students  are  U.S. 
Citizens.  A  visiting  scholar  from  China  is  also  working  on  research  related  to  this  program.  Wc 
are  now  recruiting  another  Visiting  Research  Assistant  Professor  who  can  extend  our  capabilities 
into  the  areas  of  band  structure  effects  and  quantum  transport  in  mesoscopic  systems. 


2.0  Research  Results 


2.1  Background 

Since  its  inception,  semiconductor  technology  has  been  stimulated  by  requirements  for 
electronic  systems  with  ever-increasing  capabilities  'to  process  information  faster,  more  functionally 
and  more  efficiently.  These  requirements  have  motivated  the  scaling  down  of  integrated  circuit 
(IC)  device  dimensions  into  the  submicron  (less  than  ten  thousand  angstroms)  and  ultrasubmicron 
(less  than  one  thousand  angstroms)  regions.  Today,  we  have  entered  an  era  where  nanostructure 
physics  and  fabrication  motivate  our  research  efforts  in  semiconductor  electronics  [ij.  As 
fabrication  technology  has  allowed  such  devices  to  be  realized,  many  new  and  fundamental 
questions  have  emerged  Concerning  the  underlying  physics  of  small  (atomic  level)  dimensions  in 
semiconductor  devices.  Important  issues  now  under  consideration  for  uhrasubmicron  devices 
include  nonequilibrium  transport  dealing  with  such  topics  as  quasi-ballistic  transport,  overshoot 
phenomena  and  quantum  transport.  A  great  deal  of  progress  has  been  achieved  in  our 
understanding  of  these  important  device  effects,  although  major  work  remains  to  be  done  as  our 
ability  to  fabricate  very  small  electronic  device  structures  continues  to  expand  and  mature  [2,3]. 

The  ability  to  fabricate  small  devices  has  been  continually  refined  over  the  last  fifteen  years 
through  impressive  improvements  in  materials  growth  technologies.  Molecular  beam  epitaxy 
(MBE),  metalorganic  chemical  vapor  deposition  (OMCVD),  and  atomic  layer  epitaxy  (ALE)  have 
provided  the  ability  to  fabricate  a  wide  variety  of  materials  and  hetcrostructure  combinations  with 
near  perfect  interfaces,  doping  control  and  compositional  uniformity  with  atomic  levei  dimensions. 
The  development  of  ALE  may  very  well  prove  to  be  the  ultimate  growth  technology  since  it  allows 
the  deposition  of  one  monolayer  of  device  quality  material  through  a  controllable,  self-limiting 
mechanism,  and  is  especially  useful  for  the  deposition  of  heterojunctions  [4).  Tire  ability  to  grow 
layers  with  dimensions  of  a  few  angstroms  opens  the  realm  of  quantum  transport  to  experimental 
study  and  verification.  Thus,  topics  resulting  from  size  quantization  in  condensed  matter  must  Dc 
investigated  from  theoretical  viewpoints  with  tools  which  are  either  partially  developed  or  through 
the  development  of  new  tools  which  are  not  now  available.  Quantization  effects  arising  from 
geometrical  size  constraints,  proximity  effects  resulting  from  closely  packed  arrays  of  devices,  and 
genera!  solid-state  considerations  not  heretofore  considered  questionable  (.effective  mass 
approximation,  the  role  of  interfaces  and  the  like)  must  be  addressed  from  a  fundamental  point  of 
view.  Moreover,  from  a  device  physics  viewpoint,  it  is  desirable  to  have  a  microscopic  description 
of  the  physics  of  small  dimensions  which  is  amenable  to  phenomenological  treatment,  so  that  its 
properties  can  be  meaningfully  incorporated  into  futuristic  device  concepts  and  simulation. 


Theoretical  methods  to  address  carrier  transport  have  also  progressed  rapidly  over  the  last  ten- 
fifteen  years,  in  a  manner  similar  to  research  in  semiconductor  thin  film  epitaxial  growth 
technology.  In  fact,  this  is  a  natural  progression  in  many  ways  and  is  to  be  expected.  The 
progress  achieved  in  materials  growth  of  structures  with  quantum  dimensions  dictates  that  new 
approaches  be  developed  and  refined  to  study  quantum  transport  phenomena.  However,  a 
significant  change  in  research  direction  is  now  warranted.  Past  transport  theory  and  device 
modeling  approaches  have  relied  on  particle  or  quasi-particle  approaches  where  the  electrons  are 
treated  as  rigid  mobile  entities  which  undergo  interactions  with  the  transport  medium.  The 
treatment  of  the  interaction  often  involves  wave  concepts.  However,  the  model  is  basically  a 
particle  model.  In  the  current  regime  of  quantum  transport,  we  may  no  longer  be  able  to  consider 
the  carriers  as  particles.  It  is  quite  likely  that  their  physical  behavior  will  be  governed  either 
partially  or  completely  by  wave  phenomena  (i.e.,  electromagnetic  and  quantum  mechanical 
effects). 

The  quasi-particle  methods  attempt  to  retain  as  much  of  the  classical  formalism  as  possible  in 
order  to  be  able  to  express  results  in  terms  of  parameters  which  are  of  the  greatest  experimental 
interest,  such  as  carrier  velocity  and  diffusion  constant.  This  is  a  flexible  approach.  However, 
much  care  is  required  to  ensure  that  all  important  effects  are  properly  included  because  of  the 
approximations  involved  in  the  formulations.  On  the  other  hand,  more  fundamental  quantum 
approaches,  such  as  operator-eigenfunction  methods,  adhere  closely  to  the  actual  quantum  states 
present  in  the  device  structure  when  scattering  is  not  included.  Scattering  processes  (dissipation) 
can  be  added  using  perturbation  theory  from  quantum  mechanics.  These  techniques  can  obtain  the 
greatest  sensitivity  to  the  resulting  carrier  confinement  and  the  lattice  potentials.  However,  they  are 
relatively  inflexible  in  studying  non-linear  dynamical  properties  in  the  presence  of  strong 
dissipation,  such  as  is  present  in  the  electron-phonon  interaction  at  high  electric  fields. 

Another  approach  to  quantum  transport  relies  on  the  “integral  over  paths”  method,  originally 
proposed  by  Dirac  and  formulated  by  Feynman  [5].  Path  integral  methods  rely  on  an  influence 
functional  technique  in  which  the  source  of  the  dissipation  has  been  integrated  over  all  phonon 
modes.  This  results  in  a  model  influence  functional  where  the  phonon-scattering  dissipation  can  be 
represented  as  an  interaction  with  a  collection  of  harmonic  oscillator  modes  in  which  the 
translational  invariance  of  the  carriers  is  preserved.  The  resulting  model  includes  constant  or 
oscillatory  electric  and  magnetic  fields,  carrier  screening,  scattering  and  dissipation,  carrier 
confinement,  background  temperature  and  initial  conditions  can  be  dealt  with  as  readily  as  for  a 
free  particle  f6J. 

The  work  at  N.C.  State  University  supported  by  the  ONR  under  the  current  contract  has 
progressed  over  the  past  sixteen  years  from  the  realm  of  particle  models  to  quasi-particle  models  to 
quantum  transport  models.  Wc  have  relied  primarily  on  the  Monte  Carlo  method  to  study  and 
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solve  transport  problems  in  III-V  compound  semiconductors.  Our  techniques  for  modeling 
materials  physics  and  device  phenomena  extend  to  device  dimensions  around  one  thousand 
angstroms  [7-8],  and  we  are  confident  in  these  models  for  predicting  steady  state  and  transient 
device  effects  down  to  these  dimensions.  However,  in  order  to  remain  in  the  forefront  of  transport 
physics  and  device  research,  we  must  continue  and  increase  our  progress  into  the  realm  where 
physical  effects  can  be  studied  in  device  regions  with  dimensions  less  than  one  thousand 
angstroms. 

As  one  new  direction  for  our  research,  we  have  chosen  to  pursue  the  path  integral  method. 
This  is  in  contrast  to  other  methods  now  being  studied  and  supported  by  ONR,  such  as  the 
application  of  the  density  matrix  formalism  and  Wigner  distribution  function  approaches  including 
the  use  of  moment  equations  [9].  We  have  made  substantial  progress  in  the  path  integral  approach 
during  the  past  five  years  and  we  believe  that  path  integral  methods  can  play  an  increasingly  more 
significant  role  in  developing  an  understanding  of  quantum  transport  in  devices.  Currently,  we  are 
aware  of  only  one  other  U.S.  university  program  directed  toward  applying  the  path  integral  method 
to  semiconductor  devices.  This  is  the  program  at  the  University  of  Illinois  under  the  direction  of 
Prof.  Karl  Hess.  Our  program  is  very  complimentary  to  the  one  at  Illinois,  and  we  are  in  fact 
collaborating  with  those  personnel  through  the  use  of  the  NSF  Supercomputing  Center  and  the 
National  Center  for  Computational  Electronics  there.  One  of  our  fomier  students  who  finished  his 
Ph.D.  degree  in  1990  has  joined  Prof.  Hess’  group  and  thus  we  expect  this  interaction  to  increase. 
We  propose  to  spend  a  major  effort  collaborating  with  Prof.  Hess’  group  in  applying  the  path 
integral  to  specific  device  structures  such  as  mesoscopic  devices  [10]  during  the  next  year,  thereby 
further  establishing  its  utility  as  a  method  for  modeling  quantum  transport  in  devices.  In  addition, 
we  will  continue  to  explore  the  Monte  Carlo  transport  method  for  study  of  novel  device  structures 
where  theoretical  underpinning  is  required.  One  particular  class  of  devices  which  we  will  continue 
to  study  is  based  on  hot  electron  injection  across  potential  barriers  where  the  transport  mechanisms 
which  determine  device  behavior  are  poorly  or  incompletely  understood.  Of  great  importance  here 
is  the  fact  that  these  devices  offer  the  potential  to  compare  experiment  and  theory.  Also,  these 
devices  can  be  used  to  aid  in  our  understanding  of  quasi-ballistic  transport  [3 1.  For  example,  we 
have  explored  the  concept  that  ballistic-like  behavior  can  be  predicted  even  when  scattering  effects 
are  significant  ['ll].  The  type  of  scattering  and  the  range  of  carrier  energies  are  critical  in  this 
regime.  Another  novel  class  of  structures  which  will  also  be  studied  in  detail  by  the  Monte  Carlo 
method  is  based  on  the  rapidly-emerging  pseudomorphic  or  strained-layer  devices.  These 
structures  permit  extended  compositional  ranges  and,  thus,  have  a  number  of  potential  advantages 
such  as  higher  transconductance  and  channel  carrier  density.  Recent  studies  on  the  strain-induced 
piezoelectric  fields  have  opened  yet  another  possibility  in  realizing  ultra-fast  switching  devices, 
resulting  in  an  increased  importance  on  the  study  of  carrier  transport  in  these  devices.  Othei  novel 
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devices  which  will  be  analyzed  include  delta-doped  field-effect  transistors  and  quantum  well 
structures,  based  on  studies  which  have  begun  the  past  year.  For  example,  we  have  proposed  a 
new  transistor  structure  -  the  real  space  transfer  logic  transistor  (RSTLT)  -  which  has  a  predicted 
switching  speed  of  less  than  3ps.  If  verified,  this  will  be  the  fastest  room  temperature  logic  device 
yet  reported.  Finally,  we  will  incorporate  Monte  Carlo  methods  into  the  quasi-particle  approach 
based  on  moments  of  the  Boltzmann  transport  equation  with  quantum  mechanical  corrections. 
Preliminary  results  show  the  utility  of  this  approach  for  modeling  microwave  and  millimeter  wave 
devices  which  are  important  to  the  DoD’s  MMIC  program. 

We  want  to  re-iterate  the  important  interactions  which  have  developed  between  our  research 
group  and  other  researchers  during  the  last  year.  These  interactions  increase  the  impact  of  this 
program  on  the  field  of  semiconductor  device  physics  and  allow  our  research  efforts  to  be  far  more 
productive  through  increased  intellectual  efforts  and  enhanced  facilities  and  resources.  During  the 
1990-1991  funding  periods,  we  have  intensified  our  research  collaborations  with  Prof.  K.  Hess  at 
the  University  of  Illinois,  Dr.  H.  Grubin  of  Scientific  Research  Associates,  and  Drs.  M.  Stroscio 
and  G.  Iafrate  of  the  U.S.  Army  Research  Office.  These  research  collaborations  have  resulted  in 
several  joint  publications  based  on  mutual  research  interests,  expertise  and  capabilities.  There  have 
been  several  visits  between  these  laboratories  and  the  logistics  for  increased  collaborations  are 
excellent.  Also,  we  have  very  recently  initiated  an  interaction  with  Dr.  Doran  Smith  of  ETDL,  Fort 
Monmouth  to  verify  the  operation  of  the  RSTLT. 

2.2  Summary  of  Research  Results 

2.2.1  Path-Integral  Monte  CarloiPIMC)  Research  for  Ultra-Small  Device  Aonlicarions 


This  section  will  provide  a  summary  of  work  concluded  during  the  past  year  on  the  continued 
development  of  Path-Integral  Monte  Carlo  (PIMC)  methods  for  the  study  of  the  electronic 
properties  of  ultrasmall  devices.  The  application  of  PIMC  is  based  on  the  Feynman  path-integral 
formalism  of  quantum  mechanics.  The  following  summary  represents  work  done  in  collaboration 
with  Dr.  M.A.  Stroscio  of  the  U.S.  Army  Research  Office  and  Dr.  L.F.  Register  of  the  University 
of  Illinois.  Dr.  Register  obtained  his  Ph.D.  from  N.C.  State  University  and  is  now  a  Post- 
Doctoral  Fellow  at  Illinois.  The  work  deals  with  constraints  on  the  polar-optical  phonon  influence 
functional  in  heterostructures  which  arise  form  conservation  relationships  for  certain  physical 
quantities  in  the  PIMC  formulation.  The  work  is  theoretical  and  formal;  however,  it  has  more 
practical  applications  since  it  provides  a  basis  for  estimation  of  the  influence  functional  from  limited 
knowledge  of  the  spatial  confinement  of  individual  phonon  branches.  The  results  can  help  to 
expand  our  knowledge  of  electron-phonon  interactions  in  heterostruciures. 


The  Feynman  path-integral  (FPI)  formulation  of  quantum  mechanics  allows  formal  inclusion 
of  carrier  phonon  coupling  to  all  orders  in  the  coupling  potential  via  the  influence  functional 
[12,13].  For  this  reason,  the  FPI  has  been  extended  to  the  study  of  confined  carriers  in  ultrasmall 
heterostructure  geometries  [18],  including  recent  advances  in  numerical  methods  [19-21]. 
However,  in  heterostructures,  beyond  the  confinement  of  carriers,  the  confinement  of  phonons 
should  be  considered  as  well.  Here,  this  latter  confinement  effect  on  the  influence  functionals  for 
polar  coupling  of  carriers  to  longitudinal-optical  (LO)  phonons  is  considered,  without  assuming 
any  specific  functional  form  for  the  phonon  modes.  In  Section  (a),  the  influence  functional  is 
defined  in  terms  of  the  carrier-phonon  modes.  In  Section  (b),  the  influence  functional  is  defined  in 
terms  of  the  carrier-phonon  interaction  Hamiltonian.  In  Section  (c),  it  is  shown  that  the  partial 
contributions  to  the  total  influence  functional  from  the  individual  branches  of  the  LO-phonon 
spectrum  are  related  and  constrained  by  conservation  relations.  In  Section  (c),  it  is  also  shown  that 
these  conservation  relations  lead  to  the  familiar  results  of  Feynman  [14,15]  in  the  bulk  crystal-limit, 
and  that  for  heterostructures  they  provide  a  basis  for  estimation  of  the  influence  functional  from 
limited  knowledge  of  the  spatial  confinement  of  the  individual  phonon  branches. 


(a)  Path-Integral  Formalism 

Consider  a  coupled  carrier-phonon  system  with  quantum  mechanical  Hamiltonian  of  the  form 
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where  r  is  the  carrier  coordinate,  are  the  phonon  coordinates,  fi(£>m  are  the  uncoupled  phonon 
energies,  and  <|)„/r]are  real  functions  characterizing  the  carrier  position  dependence  of  the  carrier- 
phonon  interaction.  For  this  system,  the  equilibrium  density  matrix  for  the  carrier  averaged  over 
the  phonon  coordinates  and  subject  to  a  Maxwell-Boltzmann  distribution  in  energies  at  temperature 
T  is  given  by 
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in  the  FPI  formalism  [13].  Here,  “  J...  C9,\(t)”  signifies  an  integral  over  all  paths  from  r'to  r that 
are  continuous  functions  of  time, r and  7lr(t)}  is  the  influence  functional,  accounts  for 

coupling  between  the  carriers  and  phonons  to  all  orders  in  the  coupling  potential  and  is  given  by 
f  15] 
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Equilibrium,  rather  than  time-dependent,  conditions  are  considered  here  for  compatibility  with 
existing  numerical  methods  for  performing  carrier  real  self-energy  calculations  in  confined 
geometries  [19,21  j.  However,  the  calculations  performed  here  can  be  adapted  easily  to  the  time- 
dependent  case. 


(b)  Conservation  Relations 

For  polar  coupling  of  carriers  to  phonons,  LO  or  otherwise,  the  interaction  Hamiltonian  for 
each  phonon  mode  is 
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where  £  is  the  free-space  dielectric  constant;  z  is  the  carrier  charge;  r[,  5V^and  Z[ are  the  equilibrium 

position,  mass  and  effective  charge  (the  measure  of  ionization)  of  the  /th  ion,  respectively; 

are  the  normalized  classical  lattice  vibrational  modes  of  the  crystal  lattice  as  a  discrete  function  of 

the  mode,  ion,  and  direction,  respectively;  the  subscript  j  on  a  vector  indicates  the  dot  product  of 
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that  vector  with  the  unit  vector  in  the ;th  direction;  and  a;n  and  am  are  th  raising  and  lowering 
operators,  respectively  [22].  Thus,  in  terms  of  the  vibrational  modes  <J>„,  <  dv*  coupling  functions 
<t>m(r)  of  Eqs.  (3)  and  (5)  are 
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where  the  relationship  =  V  h/2com  (a^  +  am)  has  been  used.  Of  course,  knowledge  of  (j>m(r)  is 
still  limited  by  knowledge  of  of  which  here  only  the  following  is  claimed:  for  any  crystal 
structure,  the  classical  vibrational  modes  form  a  mathematically  complete  set  of  functions  over  the 
lattice  degrees  of  freedom  and,  thus  tire  subject  to  the  orthogonality  condition 


/b,n ./  y^b/n,/  ,y  =  5/,/ ’Sj  j'  (6) 

m 

where  8  is  the  Kronecker  delta  function  [23j. 

Having  obtained  the  coupling  functions  $m(r)  in  terms  of  the  classical  vibrational  modes 
Eq.  (5),  the  same  now  can  be  done  for  the  influence  functional  ’J[r(i)l  of  Eq.  (3).  For  any 
single  branch  of  the  LO-phonon  spectrum  v  of  near  uniform  energy  htOv,  the  spatial  correlation  of 
(!)mvv(r)  encountered  on  the  riglu-hand-sidc  (RHS)  of  Eq.  (3)  becomes 
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Though  Eq.  (6)  is  exact  only  for  a  sum  over  all  modes,  for  an  entire  branch  of  the  phonon 
spectrum,  including  perhaps  noncoupling  phonon  modes  that  can  therefore  be  included  and 
weighted  as  desired,  it  is  reasonable  to  expect  the  correlation  function  Em  Om>i  j  to  remain 

short  range  in  the  displacement  n  ■  q.  Neglecting  this  displacement  relative  to  r  -  f  and  accepting 
an  uncertainty  in  ion  position  q  of  a  single  lattice  primitive  unit  cell,  transforms  Eq.  (7)  to  the  form 
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Here,  $i=  r/lr-r'l,  X  -  r'  /lr  -  r'l  and  =  rj  /I  r  -  r  1/  crion(r)  is  the  local  ion  density;  and  “<>” 
indicates  an  average  over  ions  1  within  a  lattice  primitive  unit  cell  centered  at  r,  and  lattice 
configurations  in  the  case  of  alloys.  Substituting  the  approximation  of  Eq.  (8)  into  Eq.  (3),  the 
influence  functional  takes  the  form 
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where  all  effects  of  phonon  localization  are  isolated  in  the  form  factors  rv(r/r')  [24], 

As  for  §m(r)  before,  the  result  of  Eq.  (11)  for  the  influence  functional  7lr(t)j  is  limited  by  the 
degree  of  knowledge  of  the  lattice  vibrational  modes  <J>  Fortunately,  employing  the 
orthogonality  condition  of  Eq.  (6)  exactly,  it  is  demonstrated  readily  that  the  functions  y yrjlj  (r) 
obey  the  conservation  relation 
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is  a  material  dependent  quantity  only.  For  Fv(r,r'),  via  Eq.  (9)^  this  conservation  relation 
translates  to 
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Further,  as  defined  in  Eq.  (13),  y(r)is  just  a  spatially  varying  version  of  the  usual  density  function 
encountered  in  conventional  bulk  polar  optical  scattering  rate  calculations, ;that  for  a  pure  binary 
crystal  can  be  obtained  phenomenologically  from  y1  =  co2e  (e^1  —  e’1),  where  £<*>  and  e0  are  the 
relative  high  and  low  frequency  dielectric  constants,  respectively  [25].  For  an  alloy,  a 
compositionally  weighted  averaged  of  the  pure  material  values  should  provide  a  reasonable 
approximation,  therefore,  with  these  conservation  relations,  if  not  for  the  explicitly  frequency 
dependent  terms  in  Eq.  (11),  knowledge  of  the  material  dependent  quantity  y(r)  alone  wouldkbe 
sufficient  to  completely  determine  the  total  influence  functional.  As  is,  these  conservation  relations 
still  provide  a  basis-for  estimation  of  the  influence  functional  from  limited  knowledge  of  the  spatial 
confinement  of  the  individual  phonon:  branches. 


(c)  Examples 

For  comparison  with  previous  results,  consider  the  limiting  case  of  the  influence  functional  of 
a  pure  binary  bulk  semiconductor.  For  this  case  there  is  only  one  LO-phonon  branch  and,  thus, 
the  influence  functionafis  completely  determined  by  only  that;phonon  branch’s  optical  phonon 
energy  and  the  conserved  quantity  F(r,  r')  [26].  For  a  homogeneous  material,  from  Eq.  (15),  this 
latter  quantity  is  given  by  simply  r(r,  r')  =  y(r)  =  y.  Employing  the  phenomenologically  obtained 
value  for  y  given  above,  gives  from  Eq.  (11)  exactly  the  influence  runetk/  v  •'‘blamed  by  Feynman 
for  coupling  to  polar-optical  phonons  in  bulk  semiconductors  (after  switching  to  SI  units  and 
removing  the  normalizations)  [14,  15J.  This  agreement  results  despite  the  use  of  dissimilar 
physical  models.  :Here,  discrete  dipoles  and  a  finite  Brillouin.zone  are  considered.  In  contrast, 
Feynman  employs  the  FrolichYnodel  of  the  lattice  polarization  for  carrier-polar-optical-phonon 
coupling,  that  of  approximating  the  crystal  lattice  as  a  continuum  of  dipole  charge  characterized  by 
the  density  function  y  and  the  vibrational  frequency  0)  [25],  and,  consistent  with  this  continuum 
approximation,  an'infinite  Brillouin  zone. 

As  a  second  example,  chosen  to  illustrate  a  range  of  phonon  confinement  possibilities, 
consider  an  AlGaAs-GaAs  heterostructure  of  arbitrary  geometry.  The  total  dipole  charge  densities 
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within  both  materials  are  readily  available  as  described  above.  Therefore,  given  the  conservation! 
relation  of  Eq.  (12),  only  estimates  of  the  relative  contributions  of  the  individual  phonon  branches 
as  a  function  of  position  are  required  for  evaluation  of  the  influence  functional  via  Eqs.  (9)  and 
(d  l)  [27 J.  To  this  end,  for  large  energy  separations  between  tlv-  LO-phonon  branch  in  the  GaAs 
and  the  AlAs-like  LO-phonon  branch  in  the  AlGaAs,  the  peneo  m-  >  'Vptiuof  the  AlAs-like  phonon 
modes  into  the  GaAs  should  be  small  and  in  a  first  approxir.v:  ■  can  be  neglected.  For 

small  energy  differences  between  the  LO-phonon  branch  m  *'*••  -  aAs  and  the  GaAs-like  LO- 
phonon  branch  in  the  AlGaAs,  the  branch  in  the  AlGaAs  [28  .  no  distinction- need  be  made 
between  these  two  branches,  and  thus,  the  degree  of  phonon  con'  .lent  between  these  branches, 
becomes  a  mute  point.  Therefore,  within  the  GaAs  only  .O-phonon  branch  need  be 

considered  and,  thus,  the  total  dipole  charge  density  YgoAs  *ieed  nabe  subdivided^  Further,  while 
two  phonon  branches  still  must  be  considered  within  the  AlGaAs,  the  relative  contributions  of  each 
phonon  branch  YAiAs-iikc  and  YGaAs-likc  to  the  total  dipole  charge  density  YAlGaAr.are  the  same  as. 
those  required  for  conventional  alloy  scattering  rate  calculations  [22,29-33).  If  the  energy 
separation  between  the  optical  phonon  branches  of  the  adjacent  materials  is  considered  to  be  in  an 
intermediate  range,  that  is  when  it  is  too  large  to  be  neglected  but  too  small  to  prevent  significant: 
penetration  of  the  phonon  modes  through  the  material  interfaces  on:  the  spatial  scale  of  inter"?*,  still: 
only  estimates  of  the  relative  contributions  of  each  phonon  branch  near  the  interfaces  to  the  totals 
‘/GaAs  and  YaiGoAs  arc  required,  not  absolute  nor  mode  by  mode  estimates.  Possiblefonses  for  such 
estimates  are,  for  example,  calculations  of  conf  ~cd  phonon  modes  and  penetration  depths  of 
evanescent  phonon  modes,  e.g.,  Refs.  34  -  37. 

(d)  Conclusion 

The  Feynman  influence  functional  for  polar  coupling  of  carriers  to  confined  phonon  modes  in 
heterostructures  has  been  considered.  It  was  found  that  the  partial  contributions  to  the  total 
influence  functional  from  the  individual  branches  of  the  optical  phonon  energy  spectrum  are  related 
and  constrained  by  conservation  relations,  Eqs.  (12)  -  (15).  The  derivation  of  these  conservation 
relations  required  no  specific  functional  form  for  the  phonon  modes;  rather  it  employed  only  the 
inherent  orthogonality  and  mathematical  completeness  of  the  classical  vibrational  modes  over  the 
crystal  lattice  degrees  of  freedom.  In  the  bulk  crystal  limit  these  conservation  relations  lead  to  the 
familiar  results  ofFeynman  for  heterostructures  and  these  conservation  equations  provide  a  basis 
for  estimation  of  the  influence  functional  ffom  limited  knowledge  of  the  spatial  confinements!  -he 
individualtphonOn  branches.  Further,  though  the  FPI  formalism  is  considercd  here,  it  should  be 
possible  to  apply  similar  reasoning  to  more  conventional  scattering  rate  calculations  in 
heterostructures,  as  was  done  for. optical  phonon  scaiteringtimalloys  previously  [22]-. 
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2.2.2  A  New  Non-Parabolic  Hvdrodvn;<mic  Transport  ModeUin  Semiconductors 


During  the  past  year,  we  have  made  a  major  step  towards  the  development  of  a  more  realistic 
hydrodynamic  transport  model  for  semiconductors  which  includes  non  parabolic  band  Structure 
effects.  This  model  contains  a  much  improved  physical  description  of  device  phenomena  in 
nanostructures  than  does  current  drift  diffusion  based  models  while  it  is  much  less  computer¬ 
intensive  than  Monte  Carlo  device  simulation  models.  Our  approach  has  been  accepted  for 
publication  in  the  Physical  Review  (see  Section  2.3).  This  section  will  describe  the  details  of  the 
new  hydrodynamic  model  aid  compare  the  results  wi*h  the  Monte  Carlo  method: 

The  model  is  based  upon  a  unique  derivation  of  the  moments  of  the  Boltzman'  transport 
equation  for  the  collisionless  case  (collisions  are  introduced  later).  This  derivation  implements  an 
efficient  and  compact  mathematical  formalism  appropriate  for  electrons  under  the  influence  of  high 
electric  fields  and  nonstationary  conditions.  The  theoretical  investigation  also  introduces  a  novel 
distributional  form,  with  nonparabolic  properties,  to  precisely  define  the  resulting  nonparabolic 
transport  parameters.  The  final  set  of  model  equations  are  exhibited  in  a  fashion  to  clearly  show 
the  new  factors  as  corrections  to  the  mere  familiar  hydrodynamic  model  applicable  for  the  constant 
effective  mass  case.  In  addition,  results  of  an  extensive  investigation  of  the  assumed  ansatz 
distribution  and  resulting  nonparabolic  model  parameters  using  an  elaborate  Monte  Carlo  model  are 
described.  The  Monte  Carlo  technique  was  used  to  generate  comparison  electron  distributions  and 
exact  values  for  the  nonparabolic  transport  parameters  for  stationary  and  nonstationary  electronic, 
structures.  In  all  cases,  excellent  agreement  was  found  between  the  Monte  Carlo  calculated 
parameters  and  the  derived  ronparabclic  model  terms.  The  Monte  Carlo  calculations  also  revealed, 
that  the  ansatz  distribution,  represented  a  significant  improvement  over  the  more  familiar  displaced 
Mawellian,  Therefore,  this  new  model  should  prove  very  valuable  for  studying  novel  electronic 
device  structures  operating  under  high  bias  conditions. 

(a)  General  Moment  Equations 

Consider  the  collisionless  Boltzmann  transport  equation  (BTE)  in  terms  of  the  phase-space 
variables  (r,kjt).  The  developments  presented  here  are  independent  of  the  specific  scattering 
phenomena  which  will  be  addressedtiater.  Hence,  the  BTE  is 

df  7 

j +  “■  9  * rf  + 1  *  =°  (0 

where /fr,6f)  is some  arbitrary  electron  distribution  function,  Invoking  the  effectivc-mass-theorem 
[381  (also  called  the  quasi-free-particlc  approximation  [891)  for  electron  motion  defines  the  electron, 
momentunvas  p  =  ftfoind  identifies  thc  electron  group  velocity  as 
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=  u(k)  =  jV  £'E(k).  (2) 

The  driving  force  on  the  electron  is  F  =  dp/dt  and  the  final  form  of  Eq.  (1)  is  established. 

The  overall  goal  is  to  study  moments  of  the  BTE,  Eq.  (1),  for  the  case  of  nonparabolic 
conduction  bands.  Thus,  the  first  step  is  to  develop  a  general  moment  equation.  This  is  achieved 
by  multiplying  Eq.  (1)  by  an  arbitrary  vector  (or  scalar)  operator  <b(k)  and  integrating  with  respect 
to  k  to  obtain 

|  J  ®  J  0.  (3) 

By  noting  that  <J>  is  a  function  of  k  alone  and  invoking  the  assumption  that  /  approaches  zero 
rapidly  enough  in  the  limits  of  integration  such  that  j  V  negligible  [40],  one  obtains 

the  general  moment  equation 

J;  (4)i 

in  terms  of  the  driving  force  Fj  k-space  dependent  electron  velocity  u(k)  and  an  arbitrary,  vector 
moment  operator  0(jQ.  Equation  (4)  represents  the? mathematicallform  feqoiredlto  develop  the 
nonparabolic  conservation  or  hydrodynamic  transport  equations. 

For  this  nofiparabolic  model,  the  Kane  energy  band  dispersion  relation  [38]  for  nonparabolic 
bands  v/ill  be  utilized.  This  relation  is  given  by 


(2  £2 

Zm-~'E(U(l  +  am).  (5) 

The  specific  choice  of  operators  will  determine  what  average  physical  quantity  is  conserved  in  each* 
moment  equation.  For  this  work  we  have  chosen  the  following  operators  for  use  in  the  generalized 
moment  equation,  eqn.  (4)  above. 


where 


and 


%((0  =  i 

(6) 

****•»- 

(7) 

$  2  (V  =  '£<  (k)  = “(k)  •  u(k) 

(8) 

.  I  ax/i‘6 

m(k)  =  m  %/  l  +  _m- 

m 
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‘Es(k)  =  'EM 


1  +OL'E(k)', 

a  +  2atE(H)j 


These  moment  operators  are  not  the  standard  ones  traditionally  used  in  the  development  of  the 
hydrodynamic  model.  Their  uniqueness  and  significance  are  discussed  in  detailrin  ref.  41. 

By  applying  (6),  (7),  and  (8)  in  (4),  the  nonparabolic  hydrodynamic  collisionless  transport 
equations  for  conservation  of  particles,  average  momentum  and  average  energy  in  Lagrangian  form 


dh  . 

~  =  — Vr  •  ( nv j, 
dt 


where; 


dv  5  1  __ 

—  =  - — v  •  Vrv  +  —n—  '  V  r  •  [TyJ, 
fit  r  m  nm  r  Vl 


dvJ  1 

—  =  —v  •  Vtxd+7  •  v—~  VT»(v  •  Il’-J  +  q) 
dt  n 


nv-  juffyfd3^ 


mv=  J-“  u(fQ  *  u(fOjcPfi 


«  =  -m'  f(u(V— v)  (u(U^v)fd3tj 
!T-J  =  j  mMuM  -  v)(u(k)—  v)f<T>t 

<7  =  J  Mkl  —  v)(u(k)—v)  •  (n(k)  —  v)  frt-% 


Before  this  system  of  equations  can  be  applied  to  the  study  of  electron  transport  additional 
assumptions  must  be  invoked  to  expressf-qs.  (17)  -  (-20)rtn  terms  n,  wand  w.  This  mathematical 
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closure  entails  the  use  of  specificYundamental  physical  properties  to  evaluate  theseparameters. 
Efforts  have  been  made  to  express  these  parameters  in  a  mathematicalform  which  is  as  compact  as 
possible.  It  is  clear  that  the  integral  definitions  consist  primarily  of  terms  involving  uffaand  differ 
from  the  definitions  for  the  parabolic  case  only  by  the  introduction  of  the  nonparabolic  mass  m(k) 
and  the  specific  definition  of  the  particle  velocity  u(k}.  This  property  will  be  utilized  in  the  next 
section  to  develop  approximate  expressions  so  thatithis  nonparabolic  form  of  the  hydrodynamic 
equations  may  be  used  to  study  electron  dynamics  in  realistic  electronic  structures. 


(b)  Nonparabolic  Transport  Model 

The  previous  equations  represent  a  general  nonparabolic  hydrodynamic  transport  model  for 
electron  transport  in  a  single  valley  of  GaAs.  Additional  assumptions  or  relationships  are 
necessary  to  close  the  system  of  equations  before  they  can  be  used  to  describe  transport 
phenomenon.  Specifically,  one  must  either  ignore  one  or  more  higher  order  terms  in  the 
formulation  [42],  have  some  previous  knowledge  for  the  equations  of  state  [43,44],  or  assume 
some  form  of  the  electron  distribution  function  in  momentum  space  [45],  These  are  all  acceptable 
methods  of  simplification  and  can  lead  to  the  same  unique  system  of  equations. 

The  use  of  a  predefined  distributional  form  is  often  viewed  as  being  more  restrictive  and  less 
general  than  some  of  the  previous;alternatives.  Yet,  many  derivations,  which  seek  to  develop 
models  applicableTor  more  generalized  distributions,  partially  implement  closure  relations  with 
propertiesproven  true  only  for  specific  distributional  forms.  Extending  the  electron  transport 
model  in  this  manner  has  merit  and  leads  to  results  which  areimore  physically  accurate.  However, 
mixing  assumptions  in  this  manner  can  obscure  regions-of  validity  for  the  complete  model. 
Defining  a  specific  form  of  the  distribution  function  is  the  most  complete  approach  and  easiest  to 
evaluate  from  both  a  physical  and  mathematical  point  of  view.  Therefore,  this  is  the  approach 
considered  here. 

The  distribution  most  frequently  chosen  to  describe  nonequilibrium  transport  phenomenon  is 
the  displaced  Maxwellian  function,  For  the  case  of  parabolic  conduction  bands  this  distribution 
function,  in  momentum  space,  is  [45 J 


_ uht _  f — /r  \fa--fal2] 

fj,n(k)=  u  2m'fa?t  J' 


C?l) 

L  *«»  J 

|27t7fl 

where  fa  is  the  Boltzmann  constant,  %  is  the  temperature  of  the  electron  gas  representing  the 
spread  of  the  distribution,  and  isthc  displacement  from  £=  0  of  the  Maxwellian  distribution. 
The  previous  distribution  ignores  any  affect  that  Fermi-Difac  statistics  may  have  on  the  total 
number  of  available  electrons.  Since  this  work  is  directed  toward  transport  in  conduction  bands 
and  it  is  a  simple  task  to  modify  the  normalization  coefficient  to  include  thcseTermi  effects,  they 


will  be  ignored  throughout.  This  Maxwellian  approach  has  been  applied  often  in  the  analysis  of 
such  problems;  however,  for  the  case  of  nonparabolic  bands  Monte  Carlo  calculations  indicate  that 
such  a  distribution  leads  to  a  less  than  adequate  treatment  of  the  transport  parameters  of  the 
hydrodynamic.model. 

In  this  derivation  a  unique  treatment  is  proposed  for  closing  the  general  moment  equations  for 
electron  transport  in  GaAs.  To  articulate  the  basis  of  this  approach,  consider  a  form  of  the 
displaced  Maxwellian  distribution  7dm  M  transformed  into velocity  space  such  that  7dm  (u^u  = 
fdm  (klcPk. or  specifically,  7dmM  =  Wdm  iffa’tfk))-  Here,  the  inverse  function  u'1^)  represents  the 
function  u'l(k)  =  £= flM  where  g(u)  is  the  correct  mapping  from  momentum  space  to  velocity 
space.  Also,  3dm  is  the  appropriate  Jacobian  transformation  [46]  from  momentum  space  to  velocity 
space.  Defining  this  velocity  space  distribution  in  such  a  manner  allows  parametric  definitions  of 
the  form  f=  \i(k)fdm(k)^k,  similar  to  Eqs.  (17)  through  (20),  to  be  expressed  in  velocity  space  as 
/  =  7dm  M‘Pu-  Therefore,  for  the  case  of  parabolic  bands  where  the  relationship 

between  velocity  space  and  momentum  space  is  linear,  u(£)  =  (hk)/m '  and  lJdm  •  =  this 

corresponding  displaced  Maxwellian  distribution  in  velocity  space  is 


where  v  now  represents  the  average  as  well  as  the  displacement  under  consideration  in  velocity 
space.  Both  Eq.  (21)  and  Eq.  (22)  are  symmetric  about  theirdisplacemems  within  their  respective 
spaces  and  the  average  velocity  v  is  equal  to  (hkd)M'>  therefore,  Eq.  (22)  does  represent  a 
completely  consistent  transformation  of  Eq.  (21)  for  the  case  of  parabolic  bands. 

One  would  expect  that,  under  the  conditions  of  nonparabolic  bands  where  the  effective  mass 
of  the  electron  depends  on  it’s  precise  location  in  the  conduction  band,  the  use  of  constant  m' 
would  be  insufficient.  Also,  another  dilemma  exists  in  defining  a  single  temperature  of  the  electron 
gas  since  it  has-been  shown  that  two  distinctly  different  temperatures  tensors  %  and'  %, occur  for 
the  case  of  nonparabolic  bands.  Thus,  stationary  Monte  Carlo  calcuiationsrand  physical  intuition 
suggeststhc  following  ansatz  distribution. 


7nonuar(tl)  —  n 


u,[v’r4\/Z  f—  myfv.Tjlu—vl2] 

wr-rr  2^Tw  J 


as  a  constiiuiivesrclation  to  close  the  moment  equations  for  the  case  of  nonparabolic  energy  bands. 
Here  a  single  effeciivc'icmperaiure  %,  is  still  utilized.  However, 71,,  with  a  distinct  -n subscript. 


has  been  chosen  to  replace  ^because  (3/2)%7i,  (here  *2^  is  defined  from  diagonal  elements  of  the 
energy  pressure  [TWD  approximates  the  effective  thermal  energy  well  for  the  stationary  transport 
case  at  moderate  values  of  average  energy.  In  addition,  m^y/T^bas  been  introduced  because, 
based  on  physical  considerations,  we  expect  a  nonconstant  average  effective  mass  which  should 
depend  strongly  on  the  spread  or  thermal  component  ^of  the  true  electron  distribution  and  on  the 
displacement  in  velocity  space  v. 

In  applying  Eq.  (23)  to  the  integral  definitions  (14)  -  (20),  it  is  correct  to  use  /nonpar  (u$&u  = 
f nonpar  Ikjfik,  since  /nonpar  was  postulated  in  velocity  space  from  a  distribution  which  contained  this 
transformation  property.  However,  the  issue  of  transforming  between  momentum  and  velocity 
space  for  the  case  of  nonparabolic  bands  involves  very  complicated  nonlinear  transformations. 
Thus,  it  is  instructive  to  consider  this  issue  briefly.  Assume  for  the  moment  that  a  correct  form  for 
the  distribution  in  momentum  space  is  known  and  denote  this  distribution  as  f nonpar  (%)■  To 
preserve  the  property  '/nonpar  (ukPu = so  that  transforming  integral  definitions  consists 
just  of  inverting  the  i(k}\n  :I-f&n(k)&L  a  new  nonlinear  transformation  factor,  Jnmpm  such  that 


/nonpar  (u)  ~  nonpar  I  /nonpar  (24) 

is  needed.  In  Eq.  (57),&-l-(Rl  =  £= g(u)\s  assumed  to  exist  and  to  be  a  single-valued  function  of 
u.  This  car.  be  shown  to  be  true  for  transformations  from  a  finite  region  about  k_=  0  in  momentum 
space  to  a  finite  region  about  u  =*  0  in  velocity  space:  (47J. 


For  the  case  of  nonparabolic  bands,  defined  bv  relation  (5)  and  a  resulting  nonlinear 
momentum  space  velocity  ufflgivcn  by  Eq.  (9),  the  appropriate  factor  is  1 -e-par  m where  [46J 


ch)r  <k>x 
dfc  dfy  dtU 
d\)y  d\)„  (toy 
dity  d^s 
ch)-  dv-  du- 
d*c 

Using  the  definition  of  m'# given  by  Eq.  (7)  in  Eq.  (25)  results  in 

;  j  +  c&aQ  elkfofky  citjktkz 

jnor.por  =  \^  !  I  1  *  Cf^fy  O’kjKi/z 

W  cmjL c  mm*.  i+cml 


(25) 


(26) 


cm* 


20/,- 


m  (!  +  » r  I 


(271 
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where 


The  transform  defined  by  Eq  (26)  is  extremely  complex.  Also,  it  is  obvious  that  completing 
the  remaining  determinate  operation  and  performing  the  necessary  functional  substitution  will 
introduce  cross-coupled  real  space  transformation  factors.  Therefore,  to  gain  more  insight  in  this 
investigation,  the  near  equilibrium  case  will  be  considered  where  most  of  the  electrons  are  near  the 
conduction  band  minimum.  For  this  case  the  off-diagonal  terms  may  be  neglected  and  Eq.  (26) 
can  be  approximated  to  lowest  order  in  £to  obtain 


Jnonpar 


In  order  to  proceed  further,  more  information  about  the  distribution  in  momentum  space 
appropriate  for  nonparabolic  bands  is  necessary.  The  concept  of  a  Maxwellian  function  [48], 
where  the  distribution  assumed  is  of  the  form  (n/fn)  c.\p[ — '£(k)/(k'jj%)J  (where  /„  is  the 
appropriate  normalization  factor),  will  now  be  applied  to  generate  a  test  distribution  in  momentum 
space  to  gain  additional  insight.  Doing  this  and  using  the  approximate  definition  for  electron 
energy  from  Eq.  (8)  results  in  a  simplified  i'.lustrative  distribution,  under  displacement  ktO  of 


ti 

j nonpar  Ik]  =  ~ 


— m(k. —  kj)u2  (k, — i 


Using  Eq.  (57)  as  a  test  function,  the  transformation  Eq.  (52)  yields 


n  /min-1  (k). 

7 nonpar  W  -  y  |  f_ 


kxdh2(u‘(k. —  kd)) 


ZktB'fe 


where  for  this  distribution  wl(k) represents  the  function  tr1  (k)=  k.  =  g(u) and  uv  /^represents 
some  undefined  function  of  the  form  uv  (k]  =  k, —  kd  =  8  (u'vl 

Examination  of  Eq.  (30)  reveals  that  the  exponential  argument  is  modulated  by  a  function, 
m(uv(k — kd)l  which  resembles  an  effective  mass  term  with  velocity  displacement  dependence. 

The  simplified  Jnonpar  derived  here  does  not  exhibit  the  exact  form  of  the  ansatz  distribution 
function  defined  by  Eq.  (23).  However,  the  previous  analysis  does  lend  theoretical  support  to  the 
idea  of  using  the  intuitively  derived  Jnonpar  of  Efi-  (23)  as  a  tool  to  investigate  nonequilibrium 
electron  transport  for  the  case  of  nonparabolic  conduction  bands.  Indeed,  later  it  will  be  shown 
that  this  distribution  leads  to  excellent  results  for  nonequilibrium  transport  in  nonparabolic 
conduction  bands. 

Applying  these  relations  and  through  much  mathematical  manipulation,  the  new  collisionless 
nonparabolic  hydrodynamic  transport  model  equations  are: 
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dti 

Jt 


dv 

d~t 


•  (nv) 

J  2  r  tn*  I 

lKl2;J 


dtv  _  2  m ' 

— -  =  — - v  •vrw+!F »v — •— Vr  •  | nv(u)  +  y)  (tv — ~  \v\2)\ 
at  " 


3  it 

m *  „  10a2  ,  0  ...  „„ 

| i(v,to)-  l+2u.£.(v,7it)  (u> — —  lvl2J+  3  z(v,vjjr  (n>  — ~  Ivl2/1  +  am*  \i/\2 
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The  details  of  the  derivation  are  presented  in  Ref.  41. 

The  nonparabolic  model  equations  above  have  been  written  in  a  form  so  that  for  the  parabolic 
case,  a  =  0,  it  is  easy  to  observe  that  p.  =  v  =  w  =  l  and  cp  =  0.  Then  the  parabolic  model 
conveniently  reduces  to  the  more  familiar  classical  hydrodynamic  equations  in  Lagrangian  form,  as 
given  by  eqns.  (11)  —  (20). 

(c)  Scattering  and  Transport  Model  Parameter  Evaluation 

The  nonparabolic  hydrodynamic  model  derived  in  the  last  section  ignored  any  specific  form  of 
scattering  mechanisms.  This  was  done  primarily  to  simplify  the. analysis  and  to  allow  for  a  clear 
focus  on  the  transport  parameters  which  evolve  from  the  collisionless  BTE  for  nonparabolic  bands, 
To  evaluate  the  results  of  the  previous  section  one  could  assume  a  form  for  the  scattering 
mechanism  and  then  solve  the  resulting  system  of  equations  with  appropriate  boundary  conditions 
(The  authors  have  performed  such  a  study  [50]).  Results  for  electron  current  could  then  be  used  to 
judge  the  macroscopic  feasibility  of  the  new  nonparabolic  model.  However,  the  issue  of  choosing 
accurate  forms  for  scattering  terms  is  a  complicated  one  and  the  specific  form  has  a  large  influence 
on  the  results  obtained.  Therefore,  a  different  initial  approach  will  be  followed.  The  specific 
approach  will  be  to  calculate  the  resultant  nonparabolic  parameters  m**,  [PVJ,  [P\VJ  and  q  of  the 
hydrodynamic  model  for  various  structures  using  an  elaborate  Monte  Carlo  model.  Since  the 
results  of  the  previous  section  presented  formulas  for  each  of  the  model  parameters  in  terms  of 
average  velocity  v  and  average  energy  w,  the  values  predicted  by  these  formulas  can  easily  be 
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compared  to  those  from  more  elaborate  and  detailed  Monte  Carlo  model.  Taking  this  approach  will 
give  an  excellent  gauge  of  the  accuracy  of  the  terms  of  the  model  without  the  possible  errors 
introduced  by  using  incorrect  forms  for  the  loss  mechanisms. 

The  first  test  applied  was  to  generate  the  electron  distribution  and  nonparabolic  parameters, 
directly  from  their  integral  definitions  (17)  -  (20),  for  a  time  steady-state  constant  field  structure 
over  a  variety  of  electric  field  values  using  Monte  Carlo  techniques.  For  these  stationary 
investigations,  two  independent  Monte  Carlo  models  were  used  to  analyze  electron  transport  in 
GaAs.  Both  of  these  models  utilize  single  particle  Monte  Carlo  procedures  to  improve 
computational  accuracy  and  efficiency.  The  models  incorporate  three  valley  analytical  band 
structures  and  included  all  the  significant  scattering  mechanisms  with  polar  optical  intervally 
scattering.  In  addition,  both  models  used  the  same  material  and  physical  parameters  [51]. 

Under  these  stationary  conditions,  as  can  be  observed  directly  from  Eqs.  (11)  -  (13)  or 
indirectly  in  Eqs.  (31)  -  (33),  only  the  value  of  m"  will  affect  the  final  solution  for  transport.  This 
is  true  because  all  the  other  terms  are  space  differential  in  nature.  However,  as  will  be  shown, 
these  terms  do  have  finite  values  which  vary  widely  over  the  cases  considered.  Also,  for  these 
conditions,  all  the  displacement  terms  (terms  of  form  am'  Ivl2  are  found  to  be  negligible).  Hence, 
this  is  an  excellent  test  for  the  accuracy  of  the  derived  nonparabolic  model  parameters  and  their 
dependence  on  energy  temperature  The  test  also  reveals  how  well  7W  is  determined  by  its 
dependence  on  total  kinetic  energy,  a/,  and  the  displacement  kinetic  energy  term,  (m'/Z)  \v\z.  An 
investigation  into  the  true  electron  distribution  is  also  very  important  since  Fnonpar  is  instrumental 
in  determining  the  final  formulas  for  the  nonparabolic  terms. 

The  results  of  these  stationary  investigations,  with  the  two  independent  Model  Carlo  models  in 
agreement,  are  presented  in  Figs.  1  -  4.  Fig.  1  exhibits  the  resulting  Monte  Carlo  generated 
distribution,  the  ansatz  distribution  7nonpar  and  the  displaced  Maxwellian  distribution  7dm  along 

the  direction  of  applied  field,  in  velocity  space,  for  two  values  of  electric  field.  Here,  the  Monte 
Carlo  distributions  were  obtained  by  sampling  electrons  which  resided  in  a  slice  of  velocity  volume 
from  the  three  dimensional  velocity  space.  In  the  calculations,  a  limiting  velocity  vs  is  used  to 
define  the  size  of  the  slice  in  the  directions  other  than  that  of  the  electric  field.  vs  was  chosen  small 
enough  to  clearly  and  accurately  show  the  velocity  distribution  and  large  enough  for  a  practical 
simulation  time. 

Fig.  3  clearly  indicates  that  the  corrected  electron  distribution  7nonpar  is  a  distinct  improvement 
over  displaced  Maxwellian.  Using  7dm  results  in  a  general  tendency,  for  any  select  average 

velocity  and  energy,  to  overspread  and  place  too  many  electrons  in  the  tail  of  the  distribution. 
Conversely,  the  improved  distribution  7nonpar  retains  more  electrons  closer  to  equilibrium.  This 

result  seems  correct  for  the  situation  of  an  energy  dependent  effective  particle  mass  (or 
nonparabolic  conduction  band)  and  is  verified  by  the  Monte  Carlo  generated  distributed. 
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Figure  1:  Normalized  F  valley  velocity  distribution  functions  resulting  from  the  constant  electric 
field  values  of  1.0  kVlcm and  5.0  kV/cm.  The  profiles  are  exhibited  along  the  direction  of  net 
average  velocity  (along  direction  of  applied  field).  All  distributions  shown  are  consistently 
normalized  to  contain  an  identical  number  of  electrons. 
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Figure  2:  Average  r  valley  effective  mass,  m*\  normalized  by  near  equilibrium  mass  value  of 
m*  =  0.57xlO'31kg  and  plotted  versus  average  electron  energy,  co.  These  results  were  calculated 
for  stationary  transport  conditions. 


Figure  3:  Velocity  and  energy  F  valley  pressure  tensors,  normalized  by  nks,  versus  average 
electron  energy.  The  off-diagonal  terms  are  omitted  because  they  are  negligible  (an  order.of 
magnitude  smaller)  compared  to  diagonal  terms.  These  results  were  calculated  for  stationary 
transport  conditions. 
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Figure  4:  T  valley  heat  flow  vector  along  the  direction  of  app';ed  electric  field,  normalized  by  n, 
versus  average  electron  energy.  The  other  vector  components  »u*e  omitted  because  they  are 
negligible  compared  to  field  direction  component.  These  results  were  calculated  for  stationary 
transport  conditions. 
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Figure  5:  T  valley  average  velocity  and  energy  profiles  through  a  submicron  GaAs  /Vn — N — 
N+  ballistic  diode  structure  biased  at  0.8  volts.  The  structure  consists  of  0.5  |im  long  source  and 
collector  regions  doped  1017  cm'3  with  a  0.5  pm  acdve  region  doped  10 13  cm*3.  The  diagram  also 
indicates  three  locations  ( z  =  0.53  jam, .0.85  pm  and  1.03  pm )  where  the  nonparabolic  model 
transport  parameters  were  investigated  using  the  Monte  Carlo  method. 
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Figure  6:  Normalized  average  T  valley  effective  mass,  at  three  locations  within  the  ballistic 
diode,  plotted  as  a  function  of  appled  bias. 


The  stationary  results  for  a  variety  of  electric  field  values,  corresponding  to  the  nonparabolic 
transport  parameters,  are  presented  in  Figs.  2  -  4  and  tire  plotted  Verses  the  average  electron  energy 
for  each  case.  In  each  instance,  the  nonparabolic  transport  parameters  predicted  by  the 
hydrodynamic  model  were  calculated  using  the  Monte  Carlo  values  for  average  velocity  v  and 
average  energy  w  since  they  are  fundamentally  dependent  on  these  variables. 

In  general,  the  predictions  supplied  by  the  nonparabolic  parameter  are  in  excellent  agreement 
with  the  results  calculated  by  the  Monte  Carlo  model.  Fig.  2  shows  that  the  values  for  the  average 
effective  mass  m**  predicted  by  hydrodynamic  model  are  essentially  equal  (less  than  1  percent 
difference)  to  those  calculated  by  the  Monte  Carlo  model  up  to  0.4eV.  Fig.  5,  which  gives  both 
the  momentum  pressure  [Pv);and  the  energy  pressure  [Pw],  clearly  shows  that  these  tensors  are 
approximated  well  by  a  diagonal  matrix.  Further,  the  diagonal  elements  for  both  are  in  very  good 
agreement  with  the  predictions  made  by  the  hydrodynamic  model.  Finally,  Fig.  4  shows  the 
comparison  between  the  hydrodynamic  model’s  values  for  the  heat  flow  vector  q  and  those 
calculated  by  Monte  Carlo.  These  results  also  agree  very  well,  which  leads  to  a  very  important 
point.  The  nonzero  expression  for  the  heat  flow  vector  was  derived  here  from  first  principles  for 
the  nonparabolic  conduction  band  case.  It  also  represents  the  first  nonequilibrium  analysis  for 
such  a  transport  factor  which  is  usually  based  on  thermodynamics  arguments.  However,  it  should 
also  be  noted  that  for  highly  nonstationary  conditions  this  term  may  need  to  be  determined  by 
higher  moments  of  the  BTE  before  it  is  accurate  for  transport  simulations  (52). 

While  the  previous  stationary  results  are  very  positive,  they  alone  are  not  sufficient  to  insure 
the  applicability  of  this  new  model  to  arbitrary  electron  device  structures.  Realistic  device 
structures,  operating  under  moderately  high  biases,  always  contain  regions  where  transport  occurs 
under  nonstationary  conditions.  Hence,  additional  calculations  were  performed  to  determine  the 
accuracy  of  the  nonparabolic  model  subject  to  these  conditions.  The  n+  -  n  -  n+  ballistic  diode,  a 
simple  prototypical  submicron  structure  which  has  been  extensively  studied  [53],  was  chosen  as  a 
test  device. 

A  self-consistent  ensemble  Monte  Carlo  model  was  used  to  calculate  the  nonparabolic 
parameters  for  several  strategically  selected  locations  within  this  device  as  illustrated  in  Fig.  5. 
Specifically,  locations  were  selected  which  would  exhibit  different  types  of  nonstationary 
transport.  At  location  (1)  the  electron  gas  is  relatively  cool  with  a  low  average  energy,  but,  the 
average  velocity  is  slightly  ballistic.  Converse  to  this,  at  location  (3)  the  average  energy  remains 
high  and  the  average  velocity  is  low.  Finally,  location  (2)  represents  a  case  where  both  gas 
parameters  are  elevated  and  are  experiencing  significant  gradients. 

The  results  of  these  nonstationary  studies  arc  presented  in  Figs.  6-9  where  the  various 
transport  parameters  arc  plotted  as  a  function  of  applied  bias.  Summarizing  the  results,  the 
transport  parameters  predicted  by  the  hydrodynamic  model  are  in  very  good  agreement  with  the 


20 


Average  Electron  Energy  (  eV  ) 

Figure  7:  Normalized  F  valley  velocity  pressure,  at  three  locations  within  the  ballistic  diode, 
plotted  as  a  function  of  appled  bias.  Pvlv.  represents  the  average  of  the  three  diagonal  components 
which  were  found  from  Monte  Carlo  calculations  to  be  approximately  equal.  Off-diagonal  terms 
were  negligibly  small. 


Figure  8:  Normalized  r  valley  energy  pressure,  at  three  locations  within  the  ballistic  diode, 
plotted  as  a  function  of  appled  bias.  Pw.„  represents  the  average  of  the  three  diagonal  components 
which  were  found  from  Monte  Carlo  calculations  to  be  approximately  equal.  Off-diagonal  terms 
were  negligibly  small. 
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direct  Monte  Carlo  calculations.  In  fact,  the  agreement  for  alj  parameters  are  within  10  percent  for 
the  majority  of  the  biases  studied.  In  addition,  the  difference  only  rises  slightly  higher  to 
approximately  20  percent  (for  the  heat  flow  vector)  at  low  biases  where  the  accuracy  of  the  Monte 
Carlo  statistics  are  in  question.  The  results  agreed  both  qualitatively  and  quantitatively  for  all  three 
locations  and  at  all  biases.  Hence,  all  the  results  indicate  this  new  nonparabolic  hydrodynamic 
model  can  accurately  approximate  electron  transport  in  realistic  semiconductor  devices. 

(d)  Conclusion 

This  section  has  presented  a  unique  derivation  for  a  new  nonparabolic  hydrodynamic  transport 
model.  This  new  approach  utilized  novel  moment  operators  to  reveal  simplifying  assumptions  and 
to  arrive  at  more  compact  mathematical  formulation.  This  formulation  was  then  used  to  postulate 
an  improved  ansatz  distribution  to  close  the  moment  equations  and  result  in  a  model  and 
nonparabolic  transport  parameters  which  accurately  reflects  transport  of  electrons  in  realistic 
semiconductor  conduction  bands.  When  combined  with  appropriate  dissipation  terms.  This 
simplified  model  offers  a  computationally  efficient  method  to  model  electron  conduction  in 
semiconductors.  This  type  of  model,  which  sacrifices  some  physical  detail,  leads  to  a  more 
manageable  mathematical  problem  with  the  potential  of  rapid  solution  generation  using  numerical 
methods. 

Additionally,  this  section  presents  an  extensive  Monte  Carlo  based  investigation  into  the 
applicability  of  this  new  model.  The  model  parameters  were  tested  in  both  stationary  and 
nonstationary  environments.  In  both  cases  the  results  are  excellent.  Hence,  this  model  should 
prove  instrumental  in  investigating  novel  electronic  device  structures.  The  model  offers  a  method 
to  better  approximate  electron  transport  with  the  potential  of  faster  simulation  times  as  compared  to 
the  more  computationally  intensive  Monte  Carlo  technique. 

2.2.3  Transport  in  Photoexcited  Semiconductors 

During  the  past  four  years,  we  have  studied  the  transport  of  photoexcited  carriers  in 
semiconductors.  We  have  explored  the  use  of  Monte  Carlo  simulation  as  well  as  quantum 
mechanical  wave  approaches  to  examine  scattering  effects  and  ballistic  transport.  The  specific 
device  which  we  have  chosen  to  explore  is  the  GaAs  metal-semiconductor-metal  (MSM) 
photodetector  with  a  planar  (low  electrode  capacitance)  structure  which  is  of  interest  for  discrete 
applications  and  optoelectronic  integrated  circuits  (OEICs). 

Recently,  based  on  the  results  of  our  Monte  Carlo  simulations  (sec  refs.  46,47,54  and  56  in 
Appendix  A),  Prof.  Mourou’s  research  group  at  the  University  of  Michigan  have  achieved  a 
world-record  GaAs  detector  performance.  Using  a  MSM  photodetcctor  with  0.2  |im  finger 
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spacing,  a  power  bandwidth  of  325  GHz  was  obtained.  In  last  year’s  report,  we  predicted  that 
such  high  bandwidths  were  achievable  if  parasitic  effects  could  be  minimized.  The  key  to  Prof. 
Mourou’s  results  was  the  use  of  low-temperature  MBE-grown  GaAs  to  minimize  the  effects  of 
parasitic  circuit  elements  and  hole  transport  delay  times. 

There  are  obvious  benefits  in  fabricating  MSM  photodetectors  on  GaAs.  The  high  electron 
velocity  of  GaAs  (107cm/sec)  allows  one  to  obtain  high  responsivity,  and  the  relatively  high  hole 
velocity  (3  x  106cm/sec)  makes  this  device  extremely  attractive  for  high-speed  applications. 
Advances  in  GaAs  device  fabrication  techniques  and  process  compatibility  of  the  MSM 
photodetector  with  existing  IG  processing  are  responsible  for  many  applications  of  this  device  in  its 
discrete  or  monolithic  versions.  However,  the  GaAs  photodetector  can  only  operate  at 
wavelengths  below  about  0.87  pm.  The  range  of  wavelengths  acceptable  by  this  detector  is 
unfortunately  incompatible  with  low  loss  transmission  minima  or  current  state-of-the-art  optical 
fibers.  The  low  loss  minima  of  optical  fibers  at  1.3  pm  and  1.55  jim  can  be  used  to  advantage  if 
we  investigate  semiconductors  such  as  InGaAs.  This  semiconductor  has  even  higher  electron 
velocities  than  GaAs.  However,  holes  in  InGaAs  move  much  more  slowly  than  GaAs  [54J.  This 
implies  that  the  InGaAs  detector  will  be  slower  than  GaAs  although  it  operates  at  longer 
wavelengths.  Although  some  tradeoffs  will  have  to  be  made,  there  are  interesting  areas  of 
applications  for  the  InGaAs  MSM  photodetector.  This  year,  we  have  begun  simulations  of 
GalnAs  photodetectors  for  use  at  longer  wavelengths.  The  following  section  summarizes  the 
initial  simulations  of  the  InGaAs  detector  in  an  attempt  to  define  its  performance  and  limitations. 

In  Fig.  1  we  show  a  schematic  cross-section  of  an  InGaAs  MSM  photodetector  which  we 
have  investigated.  Compositions  of  three  epitaxial  layers  are  chosen  such  that  there  is  a  perfect 
lattice  match  between  these  layers  and  the  InP  substrate.  Thus,  carrier  trapping  by  defects  created 
by  lattice  mismatch  can  be  minimized.  Since  we  already  know  that  collection  of  slow  holes 
generated  deep  beneath  the  surface  in  InGaAs  can  ultimately  determine  the  speed  of  this  device, 
the  thickness  of  the  active  InGaAs  layer  is  chosen  to  be  1  Jim.  There  is  a  minor  responsivity  loss 
associated  with  limiting  the  carrier  generation  depth  to  about  1  |im.  The  doping  densities  of  the 
active  layer  and  the  two  remaining  layers  are  dictated  by  current  progress  in  film  growth 
techniques.  The  donor  doping  densities  of  5  x  1014cnr3  for  InGaAs  and  5  x  1015cnr3  for  InAIAs 
are  very  close  to  current  background  densities  of  pure  epitaxial  layers  grown  by  MOCVD.  There 
is  no  straightforward  way  to  determine  how  thick  the  InAIAs  barrier  enhancing  layer  should  be.  A 
semiconductor  with  a  larger  bandgap  must  be  introduced  between  the  InGaAs  and  the  Schottky 
contact  to  lower  the  dark  current.  The  thicker  this  barrier,  the  smaller  will  be  the  probability  for 
carrier  tunneling  through  the  barrier.  Hence,  the  dark  current  will  be  smaller.  However  it  is  also 
true  that  a  thicker  barrier  would  tend  to  keep  photogeneratcd  carriers  inside  the  InGaAs  for  longer 
periods  of  time.  Based  on  available  experimental  designs  of  InGaAs  MSM  photodeicctors  f55|  we 
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Figure  1:  Cross-sectional  view  of  IriGaAs  MSM  photodetector.  The  toana  bottom  InAlAs  layers 
are  used  to  confine  carrier  generation  to  1  pm  thick  InGaAs.  The  broken  line  represents  the 
simulated  are  of  the  device. 
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Figure  ">.  Two-dimensional  initial  pattern  of  generated  electrons  and  holes.  There  are 
2000  generated  electron-hole  pairs. 


chose  500 A  for  the  thickness  of  the  InAIAs  barrier  enhancing  layer.  A  separate  study  will  be 
necessary  to  determine  the  optimal  barrier  thickness.  A  sub-half  micron  distance  betweerifingers 
is  chosen  since  there  are  no  known  simulation  results  for  such  devices.  All  present  photodetectors 
are  usually  fabricated  with  1-3  |im  distance  between  fingers.  The  simulation  technique  employed 
in  the  present  studies  does  not  permit  the  simulation  of  long  devices  due  to  a  prohibitively  long 
^snputation  time.  Therefore,  the  presented  results  are  fora  0.2  pm  detector. 

The  initial  electron-hole  pair  distribution  in  InGaAs  is  determined  by  the  finger  spacing  and 
wavelength  of  light.  Since  the  spacing  between  fingers  is  smaller  than  the  wavelength  of  light, 
there  will  be  a  diffraction  pattern  formed  in  the  interfinger  area.  We  approximate  this  pattern  with  a 
Gaussian  exponentially  decaying  with  depth  as  shown  in  Fig.  2.  The  exact  pattern  could  be 
Obtained  from  Maxwell’s  equations.  The  absorption  coefficient  for  InGaAs  used  in  generating  the 
initial  carrier  pattern  is  1.25  x  104cm_1  [56].  The  number  of  electron-hole  pairs  used  in  the 
simulation  is  2000.  We  have  performed  simulations  with  various  initial  energies  of  electrons  and 
holes  or,  alternatively  various  wavelengths  of  incident  light. 

The  intrinsic  time  response  shown  in  Fig.  3  is  obtained  from  simulations  with  a  wavelength  of 
1.55  mm.  Thus,  the  electrons  are  gencrated  with  an  initial  energy  slightlyriarger  than  the  InGaAs 
bandgap.  This  curve  displays  similar  features  as  the  time  responses  for  GaAs  detectors,  except  for 
the  change  in  time  scale.  Since  the  initialresults  revealed. that  it  is  necessary  to  increase  the  bias 
voltage  to  increase  the  detector  speed,  this  simulation  was  performed  for  a  bias  voltage  of  1.4  V. 
The  InAIAs  layer  has  been  removed  for  these  simulations.  We  see  that  this  device  is  considerably 
slower  than  the  GaAs  device.  As  we  have  done  for  the  GaAs  photodetectors,  discrete  and 
monolithic  versions  of  the  detector  are  introduced.  By  comparing.the  time  responses  of  these  two 
versions  of  the  photodetector,  we  notice  only  minor  differences  between  them.  This  behavior  is 
caused  by  a  much  slower  intrinsic  pulse  for  InGaAs  than  for  GaAs.  We  deal  with  a  parasitic 
circuit  suitable  for  high  frequencies,  or  suitable  for  generating  short  pulses,  excited  by  a  relatively 
slow  current  source.  As  a  consequence,  we  cannot  take  advantage  of  monolithic  integration  of  this 
device  since  reduction  in  values  of  the  parasitics  does  not  have  much  influence  on  the  response. 
We  see  from  Fig.  3  that  it  takes  much  longer  than  25  ps  for  the  output  voltage  to  settle.  The  long 
tail  associated  with  holes  causes  the  bandwidth  of  the  detector  to  decrease. 

In  Fig.  4  we  show  the  output  power  spectra  of  an  InGaAs  MSM  photodetector.  One  common 
feature  of  ail  curves  is  that  for  frequencies  up  to  about  100  GHz.  there  is  little  difference  in  output 
power  of  intrinsic,  d:screte,  and  monolithic  detectors.  Consequently  the  3  dB  bandwidth  of  2-1.8 
GHz  for  the  intrinsic  detector  remains  almost  unchanged  for  the  discrete  and  monolithic  versions. 
We  can  make  the  following  summary  based  on  these  initial  results  for  the  InGaAs  MSM 
photodctector.  The  device  is  useful  for  detection  of  optical  pulses  at  1.55  pm.  1.  is  necessary  to 
optimize  the  thickness  of  the  barrier  enhancing  layer  (InAIAs)  in  order  to  lower  the  dark  current 
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Figure  3  Simulated  intrinsic  and:extrinsic:  response;  for  an  InGaAs  photodetector. 
Much  slower  holes  in  InGaAs  than  in  GaAs  cause  this  device  to  have  a  long  hole  t2il.  This 
also  results  in  very  similar  responses  for  the  discrete  and  monolithic  versions  of  the  device. 


Frequency  (GHz) 

Tigurc  4  Power  spectra  of  InGaAs  MSM  photodetector.  The  3  cB  frequency  «s  21.S 
GHz  for  ait  considered  versions. 


without  significant  deterioration  of  electron  collection.  Although  the  small  distance  between 
fingers  suggests  possible  use  of  small  bias  voltage,  there  is  a  speed  deterioration  associated  with 
small  bias.  In  order  to  remove  slow  holes  and  increase  the  device  speed  we  must  increase  the  bias 
voltage,  which  unfortunately,  increases  the  dark  current.  The  interdigitated  pattern  shows  similar 
responses  for  the  discrete  and  monolithic  detectors. 

2.2.4  Study  of  Delta-doped  FET  Structures 

(a)  Monte  Carlo  Simulation  of  Delta-Doped  HEMT  Structures 

As  a  promising  candidate  for  ultradiigh  speed  and  microwave  applications,  delta-doped  FET 
structures  are  attracting  considerable  attention.  By  applying  the  delta-doping  technique,  devices 
(such  as  MESFETs,  AlGaAs/GaAs  based  HEMTs,  pseudomorphic  HEMTs,  and  HBTs)  can 
achieve  improved  performances.  These  include  enhanced  threshold  voltage  control,  reduced 
trapping  effects,  better  breakdown  characteristics,  and  (most  importantly)  increased  conducting 
channel  electron  density  and  more  efficient  utilization  of  introduced  dopants.  The  resulting 
expc  'mental  devices  show  higher  current  drive  capability,  higher  transconductance  and  higher 
frequency  performance.  In  spite  of  continuing  experimental  interests  in  delta-doped  structures, 
theoretical  studies  on  these  devices  are  rare  and  mainly  first  order  (either  analytical  model  or  by 
drift-diffusion  approximation).  Theoretical  study  of  electron  transport  properties  in  such  structures 
is  necessary  in  order  to  understand  the  mechanism  of  device  operation,  to  further  improve 
(optimize)  device  performance,  and  to  exploit  potential  new  device  structures.  During  the  past 
year,  we  have  performed  a  systematic  study  of  the  delta-doped  HEMT  structures.  In  particular,  the 
following  issues  are  theoretically  investigated  for  submicron  AlGaAs/GaAs  delta-doped  HEMTs 
by  ensemble  Monte  Carlo  simulation: 

•  Comparative  study  of  delta-doped  and  uniformly  doped  HEMTs  —  First,  we  compared 
the  performances  of  experimentally  available  (from  literature)  submicron  HEMTs  and  a 
comparable  uniformly-doped  HEMT  and  investigated  mechanisms  for  the  enhanced 
perform0  ,ee  of  delta-doped  HEMTs.  To  correctly  simulate  the  operations  of  practical 
devices,  our  Monte  Carlo  model  includes  the  real-space  transfer  effect,  the  carrier 
degeneracy  effect,  and  influences  of  DX  centers  and  surface  states  in  addition  to  all  the 
scattering  mechanisms  (such  as  polar  optical,  intervalley,  ionized  impurity,  electron- 
electron  scatterings).  A  two-dimensional,  self  consistent  Poisson  er  /  tlon  solver  is 
implemented  in  the  program  for  accurately  updating  the  electric  field  in  the  device.  Our 
results  reveal  that  the  improvement  in  electron  confinement  in  conducting  channel  (as 
shown  in  Fig.  1)  and  reduced  peak  channel  electric  field  are  the  main  reasons  that  product. 
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Figure  1:  Cross-sectional  view  of  conduction  band  profiles  at  Vgs  =  0.4  V  and  V(jS  =  2.0V.  The 
cross  section  is  taken  -perpendicular  to  the  middle  of  the  0.5  (im-long  gate.  The  zero  of  the 
coordinate  (i.e.,  0  A)  is  chosen  at  the  gate  contact.  An  A1  mole  fraction  of  0.3  and  a  Schottky 
barrier  height  of  0.8  e\  are  assumed  for  both  the  delta-doped  HEMT  (D-HEMT)  and  the  uniformly 
doped  HEMT  (U-HEMT). 
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Figure  2:  Transconductance  as  a  function  of  gate-to-source  bias  VgS  at  Vds  =  2.0  V  in  a  0.5-M.m 
gate-length  delta-doped  HEMT:  (a)  For  cases  B,C,  and  E,  the  DX  center  concentration  ratio  varies 
from  0  (B)  to  0.2  (C)  to  1.0  (E)  while  the  surface  state  density  is  Fixed  to  8.0  x  10l°  cm*2,  (b)  For 
cases  F,  C,  and  G,  the  surface  state  density  varies  from  0  (F)  to  8.0  x  1010  cm-2  (C)  to 
2.0  x  1012  cm'2  while  the  DX  center  concentration  ratio  is  Fixed  to  0.2. 


reduced  parallel  conduction,  higher  electron  drift  velocity  and  increased  transconductance 
for  the  delta-doped  HEMTs,  especially  at  high  (positive)  gate  bias.  We  show  in  detail 
how  band  bending  in  the  two  structures  changed  the  electron  occupations  in  different 
layers  within  the  devices  (that  is,  top  AlGaAs  layer,  channel  layer,  and  substrate  layer). 
Our  simulation  results  also  indicate  that  even  though  the  delta-doped  HEMT  has  less  total 
AlGaAs  doping  than  uniformly  doped  HEMT  for  the  same  threshold  voltage  and 
capacitance,  performance  improvement  for  delta-doped  HEMTs  can  be  realized  due  to  the 
doping  profile  design.  The  advantage  of  delta-doped  HEMTs  over  its  uniformly  doped 
counterparts  in  terms  of  high  channel  electron  density,  reduced  parallel  conduction  in 
AlGaAs  layer,  higher  current  and  transconductance  are  clearly  demonstrated  in  this  study. 

Investigation  of  DX  centers  and  surface  states  on  delta-doped  HEMTs  —  The  roles  of 
DX  centers  and  surface  states  associated  with  doped  AlGaAs  layer  are  studied  in  this 
phase  of  our  work.  DX  centers  and  surface  states  are  believed  to  be  the  main 
mechanisms  that  degrade  performances  of  HEMT  structures.  We  studied  delta-doped 
HEMTs  under  different  defect  situations  (no  defects,  DX  centers  only,  surface  states 
only,  and  both  DX  center  and  surface  states)  under  different  bias  conditions.  We 
demonstrated  that  both  DX  centers  and  surface  states  indeed  have  negative  effects  on 
delta-doped  HEMTs  and  provided  physical  explanations  for  the  performance 
degradations.  The  analysis  of  bias  dependence  on  the  influences  of  the  two  defect  levels 
indicates  that  the  degrading  effects  are  more  significant  with  positive  gate-to-source 
voltage.  We  found  that  the  reasons  for  the  degradation  under  the  presence  of  DX  centers 
and  surface  states  are  through  reduced  channel  electron  concentration,  increased  channel 
electron  intervalley  scatting,  and  enhanced  real-space  transfer.  For  the  individual  role  of 
each  defect  state,  we  demonstrated  that  DX  centers  have  more  detrimental  influence  on 
drain  current  and  transconductance  than  the  surface  states  for  delta-doped  HEMTs  (as 
shown  in  Fig.  2). 

Simulation  of  doping  profile  variation  on  delta-doped  HEMT  performance  —  Since  the 
device  improvements  for  delta-doped  HEMTs  result  from  the  delta-doping,  variation  in 
the  doping  profile  will  play  an  important  role  on  device  performance  optimization.  We 
focused  our  study  on  the  doping  profile  influence  in  this  work.  Various  doping  situations 
are  studied.  Specifically,  we  varied:  a)  the  (AlGaAs)  peak  doping  concentration  and 
doping  layer  thickness  (single  delta-doped  HEMTs)  under  constant  threshold  voltage 
conditions;  and  b)  the  spacing  between  two  delta-doped  layers  for  double  delta-doped 
HEMT  structures  with  identical  total  doping  in  the  AlGaAs  layer.  Although  the  studied 
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double  delta-doped  HEMTs  are  experimentally  unavailable,  theoretical  simulation  of  this 
device  is  both  necessary  and  important  to  predict  their  performances  and  potential 
problems.  Theoretical  study  of  this  new  structure  can  provide  valuable  guidance  for 
device  design,  fabrication,  and  characterization.  In  this  work,  device  performances  in 
terms  of  threshold  voltage  sensitivity,  current  drive,  transconductance,  cut-off  frequency, 
and  characteristic  switching  time  are  examined.  Our  simulation  results  show  that:  I)  For 
single  delta-doped  HEMTs  with  identical  threshold  voltages,  the  thinner  the  doped  layer 
(with  more  heavily  doped  concentration),  the  better  the  overall  device  performance.  2) 
For  double  delta-doped  HEMTs  with  an  identical  total  doping  in  the  AlGaAs  layer,  an 
increase  in  the  spacing  between  two  delta-doped  layers  leads  to  improved  threshold 
voltage  design  but  causes  a  slight  degradation  in  transconductance,  cut-off  frequency,  and 
switching  time.  As  the  gate  bias  increases,  this  degradation  becomes  less  pronounced 
due  to  the  higher  “onset”  of  parallel  conduction  (for  the  structures  with  a  larger  spacing 
between  the  delta-doped  layers). 

(b)  Analysis  of  Submicron  Delta-Doped  GaAs  MESFETs 

Efforts  to  achieve  high  speed  and  high  frequency  has  led  to  continued  efforts  to  develop  novel 
device  structures.  One  of  the  most  promising  structures  is  the  MESFET  which  employs  the  delta¬ 
doping  technique.  Delta-doped  GaAs  MESFETs  based  on  MBE  and  MOCVD  materials  have  been 
investigated  experimentally  by  several  research  groups.  The  results,  while  promising,  show  a 
wide  diversity.  The  reported  maximum  transconductance  for  0.5  pm  channel-length  delta-doped 
MESFETs,  for  example,  range  from  75  mS/mm  to  400  mS/mm  (as  provided  by  both  simulation 
and  experimental  device  measurement  in  this  work).  Studies  show  that  maximum  intrinsic 
transconductance  of  the  delta-doped  MESFET  is  higher  than  that  of  comparable  HEMT  in  the 
channel  length  range  of  0.1  to  2.0  pm.  Differences  between  experimental  results  and  theoretical 
predictions  may  stem  from  material  properties  and  processing  conditions.  On  the  other  hand, 
design  and  optimization  of  such  devices  play  an  important  role  in  the  improvement  of  device 
performance.  It  is  of  practical  importance  to  identify  the  key  parameters  for  optimized  device 
performances  and  to  establish  trade-offs  needed  for  specific  applications.  It  is  obvious  that  a 
comprehensive  study  is  required  to  investigate  the  influences  of  delta-doped  MESFET  parameters 
on  performance  and  provide  useful  design  guide. 

In  this  work,  we  used  a  two-dimensional,  drift-diffusion  model  to  simulate  the  electrical 
characteristics  of  submicron  delta-doped  GaAs  MESFETs.  Poisson  and  electron  current  continuity 
equations  were  solved  self-consistentiy  and  an  analytical  velocity-field  relation  which  accounts  for 
the  velocity  overshoot  was  employed.  The  results  were  compared  with  experimental 
measurements  and  showed  very  good  agreement.  Material  parameters  are  based  on  the  results 
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from  delta-doped  MBE  materials  provided  by  Dr.  Bedair’s  group  in  the  NCSU  Electrical  and 
Computer  Engineering  department.  The  device  results  were  from  DC  measurement  data  provided 
by  Texas  Instruments  based  on  the  above  materials.  In  this  work,  we  studied  the  following  issues: 

•  Delta-doping  profile  —  We  used  a  Gaussian  distribution  to  approximate  the  actual  doping 
concentration  profile.  Simulation  of  devices  with  0!5  (im  channel  length  gives  very  good 
agreement  with  fabricated  experimental  devices.  For  delta-doping  profile  simulation,  we 
first  simulated  delta-doped  MESFETs  with  different  peak  delta-doping  concentrations  and 
half  width  at  half-maximum  (HWHM)  under  constant  threshold  voltage;  then  the 
influence  of  the  gate-to-channel  distance  was  studied.  We  show  that  steeper  doping 
profiles  (i.e.,  high  peak  concentration  and  small  HWHM  of  delta-doping)  and  a  short 
gate-channel  distance  are  preferred  for  improvement  in  transconductance.  Current  drive 
capability  is  increased  for  devices  with  steeper  doping  profile,  but  is  decreased  for 
devices  with  short  gate-to-channel  distance  (due  to  the  upward  shift  of  threshold  voltage 
as  indicated  by  experimental  results). 

•  Top  layer  background  doping  —  We  demonstrate  in  this  work  that  a  trade-off  should  be 
made  in  the  design  of  the  top  layer  background  doping.  Calculations  show  that 
significant  increase  in  current  drive  (~  200%)  can  be  obtained  at  a  low  cost  of  slightly 
reduced  transconductance  (~  10%)  for  devices  using  a  relatively  high  top  layer  doping 
concentration  (up  to  20%  of  the  peak  delta-doping).  The  threshold  voltage  shift  with 
respect  to  gate-to-channel  distance  for  MESFETs  with  higher  top  layer  background 
doping  is  considerably  higher  in  absolute  value  but  considerably  lower  in  relative  change 
when  compared  to  a  ‘true’  delta-doped  device  (with  unintentional  top  layer  background 
doping). 

•  Device  lateral  feature  size  —  We  performed  scaling  simulations  for  device  lateral  feature 
size.  The  result  show  that  reducing  the  lateral  features  sizes  (both  channel  length  and 
source-to-drain  spacing)  brings  improved  transconductance,  current  drive  capability,  and 
cut-off  frequency.  For  devices  with  identical  channel  length,  a  self-aligned  structure  is 
preferred  in  order  to  obtain  minimum  source-to-drain  spacing.  The  improvement 
obtained  by  reducing  the  gate  length  is  shown  to  be  more  profound  than  that  achieved  by 
reducing  the  source-to-drain  spacing.  Simulation  results  show  that  the  main  reason  for 
the  improvement  is  realized  from  an  increase  in  average  electron  drift  velocity  along  the 
device  channel  with  reduced  feature  sizes. 


2.2.5  Study  of  a  Novel  Heteroiunction  Real-Space  Transfer  Logic  Transistor  (RSTLT) 

In  this  section,  we  report  very  preliminary  results  for  a  new  heterojunction  logic  device.  This 
device  is  called  the  real-space  transfer  logic  transistor  (RSTLT).  These  results  will  be  presented  at 
the  International  Electron  Devices  Meeting  (EEDM)  which  will  be  held  in  Washington,  DC,  during 
December  8-11,  1991.  The  device  structure  consists  of  a  GaAs  channel  which  is  sandwiched 
between  two  AlGaAs  barrier  layers.  It  has  four-terminals  for  operation  -  a  source,  drain,  and  two 
collectors.  The  drain  terminal  is  the  signal  input  and  the  logic  operation  is  provided  at  the  two 
collectors.  The  basic  structure  realizes  NOT  and  EQUIVALENT  logic  operations,  and 
NAND/NOR  operation  can  be  achieved  with  two  of  the  devices.  A  conservative  estimate  of  the 
switching  speed  is  approximately  3  psec  total  propogation  delay.  This  estimate  is  based  on  self- 
consistent  two-dimensional  (ensemble)  transient  Monte  Carlo  simulations.  The  operation  of  the 
device  is  based  on  the  principle  that  the  spatial  location  for  real-space  transfer  can  be  controlled  by 
an  input  signal  which  steers  current  flow  under  one  or  both  of  the  collector  terminals.  Thus,  to 
perform  logic  operations,  the  carriers  do  not  have  to  traverse  the  comple.e  length  of  the  channel. 
The  principal  speed  limitation  is  determined  by  the  rate  of  real-space  transfer  in  the  device. 

The  principles  of  operation  and  further  discussion  of  device  properties  will  be  presented  in  a 
forthcoming  paper. 

2.3  Publications  and  Presentations 

The  funds  for  this  program  did  not  arrive  at  the  University  until  March  1991.  Thus,  we  have 
been  working  on  the  current  contract  for  slightly  more  than  seven  months.  However,  during  the 
last  year  we  have  made -six  oral  presentations  at  national  and  international  conferences.  Also, 
twelve  written  manuscripts  have  been  published  in  the  refereed  literature.  In  addition,  six  other 
manuscripts  have  been  accepted  for  publication  and  five  manuscripts  have  been  submitted  or  are 
currently  in  preparation.  The  following  paragraphs  summarize  the  presentations  and  publications 
made  under  this  program  during  the  last  year,  and  include  material  not  previously  reported  to  ONR 
in  the  last  renewal  proposal  as  well  as  material  which  was  accepted  for  publication  in  1990  and  has 
since  been  published  in  1991. 

A.  Conference  Presentations  and  Seminars 

K.W.  Kim,  M.A.  Stroscio  and  J.C.  Hall,  “Frequencies  of  Confined  Longitudinal-Optical  Phonon 
Modes  in  Short-Period  Strained  Semiconductor  Superlaitices,”  presented  at  the  SPIE  International 
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Symposium  on  Optical  and  Optoelectronic  Applied  Science  and  Engineering,  San  Diego,  CA,  July 
1990. 

W.C.  Koscielniak,  M.A.  Littlejohn  and  J.L.  Pelouard,  “Physical  Speed  Limits  of  Metal- 
Semiconductor-Metal  Photodetectors,”  presented  at  the  17th  International  Symposium  on  Gallium 
Arsenide  and  Related  Compounds,  Jersey,  Channel  Islands,  24-27  September  1990. 

M.A.  Stroscio,  K.W.  Kim  and  M.A.  Littlejohn,  “Theory  of  Optical-Phonon  Interactions  in  a 
Rectangular  Quantum  Wire,”  presented  at  the  SPIE  International  Conference  on  Physical  Concepts 
of  Materials  for  Novel  Optoelectronic  Device  Applications  II:  Device  Physics  and  Applications, 
Aachen,  Federal  Republic  of  Germany,  October  1990. 

M.A.  Stroscio,  G.J.  Iafrate,  K.W.  Kim,  M.A.  Littlejohn,  H.  Gronkin  and  G.  Maracus. 
“Transition  from  LO-Phonon  Scattering  to  SO-Phonon  Scattering  in  Polar  Semiconductor 
Superlattices,”  presented  at  the  Workshop  on  Compound  Semiconductor  Microwave  Materials  and 
Devices,  Fort  Lauderdale,  FL,  February  1991. 

M.A.  Stroscio,  K.W.  Kim  and  Amit  Bhatt,  “Electron-Optical-Phonon  Scattering  Rates  in 
Rectangular  Quantum  Wire,”  presented  at  the  Annual  March  Meeting  of  the  American  Physical 
Society,  St.  Louis,  MO,  March  1991. 

D.L.  Woolard,  R.J.  Trew  and  M.A.  Littlejohn,  “Effect  of  Contacts  on  Transit  Time  and  Velocity 
Overshoot  in  GaAs  n+-n-n+  Structures,”  presented  at  the  Workshop  on  Compound  Semiconductor 
Microwave  Materials  and  Devices,  Fort  Lauderdale,  FL,  February  1991. 

(B)  Refereed  Publications,  Papers  Accepted  for  Publication  and  Papers  Submitted  or  in 
Preparation 

K.W.  Kim,  M.A.  Stroscio,  and  J.C.  Hall,  “Frequencies  of  Confined  Longitudinal-Optical  Phonon 
Modes  in  Short-Period  Strained  Semiconductor  Superlattices,”  Proceedings  of  the  SPIE 
International  Symposium  on  Optical  &  Optoelectronic  Applied  Science  and  Engineering,  1336, 43 
(1990). 

M.A.  Stroscio,  K.W.  Kim,  and  J.C.  Hall,  “Variation  in  Frequencies  of  Confined  Longitudinal- 
Optical  Phonon  Modes  due  to  Changes  in  the  Effective  Force  Constants  Near  Heterojunction 
Interfaces,”  Superlattices  and  Microstructures,  7, 115  (1990). 

M.A.  Littlejohn,  J.L.  Pelouard,  W.C.  Koscielniak  and  D.L.  Woolard,  “Device  Simulation 
Augmented  by  the  Monte  Carlo  Method,”  Computational  Electronics:  Semiconductor  Transport 
and  Device  Simulation.  K.  Hess,  J.P.  LeBurton  and  U.  Ravaioli,  eds.,  Kluwer  Academic 
Publishers,  pp.  69-75,  1991. 

D.L.  Woolard,  M.A.  Stroscio,  M.A.  Littlejohn,  R.J.  Trew  and  H.L.  Grubin,  “A  New  Non¬ 
parabolic  Hydrodynamic  Transport  Model  with  Quantum  Corrections Computational  Electronics: 
Semiconductor  Transport  and  Device  Simulation.  K.  Hess,  J.P.  LeBurton  and  U.  Ravaioli,  eds., 
Kluwer  Academic  Publishers,  pp.  59-62,  1991. 

W.C.  Koscielniak,  M.A.  Littlejohn  and  J.L.  Pelouard,  “Physical  Speed  Limits  of  Metal- 
Semiconductor-Metal  Photodetectors,”  Gallium  Arsenide  and  Related  Compounds  1990.  Institute 
of  Physics  Conference  Series  No.  1 12,  K.E.  Singer,  editor,  Institute  of  Physics,  New  York,  1 12, 
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ABSTRACT 


The  theories  describing  the  dispersion  of  longitudinal-optical  (LO)  phonon 
modes  and  electron — LO-phonon  interactions  are  generalized  to  include  the  effects 
of  both  strain  and  confinement  in  semiconductor  superlattices  and  semiconductor 
microstructures.  In  particular,  the  effects  of  strain  and  confinement  on  LO 
phonon  frequencies  are  analyzed  for  short-period  strained-layer  superlattices  as 
well  as  for  semiconductor  microstructures  where  effective  force  constants  are 
altered  near  heterojunction  interfaces. 

1.  INTRODUCTION 


Confined  phonons  in  semiconductor  heterostructure3  have  recently  been  the 
subject  of  both  theoretical  and  experimental  investigations.  These  efforts  are 
due,  in  part,  to  the  fact  that  scattering  by  the  longitudinal-optical  (LO)  phonon 
modes  is  an  important  energy  loss  mechanism  for  a  wide  variety  of  III-V  semicon¬ 
ductor  devices.  Raman  scattering  measurements  have  been  used  to  study  the 
properties  of  phonon  modes  in  quantum  wells  and  super  lattices .  ~  In  particular, 
binary  GaAs/AlAs  superlattices  have  been  experimentally  investigated  mese 
extensively  due  to  the  simplicity  in  analysis  (compared  to  tenary  GaAs/AlGaAs 
superlattices)  and  the  material  familiarity.  It  is  well  known  that  phonon 
modes  in  the  acoustic  branch  are  propagative  across  the  heterojunctions  and,  as 
a  result,  are  ’’folded"  into  the  reduced  Brillouin  zone,  while  the  modes  in  the 
optical  branch  are  "localized"  in  each  layer  leading  to  descriptions  of  confined 
modes  or  slab  modes.  Furthermore,  the  localization  of  ohonon  modes  near  the. 
heterojunctions  leads  to  the  well-known  interface  modes.  In  addition  It  is 
suggested  chat  the  Raman  frequency  shifts  for  confined  LO-phonon  modes  can  be 
used  as  a  means  to  determine  the  bulk  phonon  dispersion  relation  (i.e.,  the 
equivalent  wave-vector  model)4  and  as  a  probe  to  characterize  the  quality  of 
interfaces.  Frequency  shifts  due  to  confinement  have  been  anslyzed  theo¬ 
retically  as  well;  the  results  obtained  by  a  simple  linear-chain  model  show  a 
remarkable  agreement  with  experimental  data  except  In  monolayer  superlattices. 
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A  linear-chain  model  with  variable  force  constants  for  the  couplings  at  heterojunc¬ 
tions  is  used  to  calculate  the  spectrum  of  confined  longitudinal-optical  (LO)  phonon  fre¬ 
quencies  in  a  short-period  strained-layer  superlattice.  Even  for  the  case  of  a  superlattice 
with  only  two  atomic  monolayers  in  each  superlartice  layer,  it  is  demonstrated  that  the 
frequencies  of  confined  LO-phonon  modes  arc  only  weakly  dependent  on  variations  in 
interfacial  force  constants. 


1.  Introduction 

Phonons  in  semiconductor  heterostnictures  have 
received  much  attention.  These  efforts  are  due.  in  pan.  to 
the  fact  that  scattering  by  the  longitudinal-optical  (LO) 
phonun  modes  is  an  important  energy  loss  mechanism  for 
electrons'  in  a  wide  variety  of  IH-V  semiconductor  dev¬ 
ices.1  In  addition,  Raman  scattering  measurements  have 
been  applied  to  study  the  properties  of  phonon  modes  in 
quantum  wells  and  superlattices.2-8  Recently,  Raman 
spectra  have  been  used  to  characterize  monolayer  varia¬ 
tions  in  the  thicknesses  of  quantum  wells,8  In  the  analysis 
and  interpretation  of  Raman  data  for  confined  LO-phonon 
modes,  it  is  convenient  to  model  the  confined  LO  modes 
as  modes  in  a  uniform  isolated  dielectric  slab.5-11  Such 
an  approximate  model  is  based  on  the  assumption  that 
atomic  force  constants  at  hetcrojunction  interfaces  are 
identical  with  those  of  the  bulk  or  of  unifurm  pscu- 
domorpltic  layers  In  this  paper,  the  effect  of  varying 
force  constants  at  the  hetcrojunction  interfaces  of  a 
strained-layer,  short-period  .supcrluiiice  is  analyzed  by 
using  a  linear-chain  model  to  ricicnuiue  the  corresponding 
changes  in  the  frequencies  of  the  confined  LO  phonon 
modes'  In  the  studies'  reported  herein,  a  supcrluiiice  with 
only  two  atomic,  monolayers  in  each  supcrluiiice  layer  is 
mudclcd  in  order  lo  deieiminc  the  iiupoiiaiiLc  of  pcnuibu 
lions  in  heicropinciion  interlace  force  constants  fin  a  sys¬ 


tem  where  such  effects  should  be  large.  In  particular,  fre¬ 
quencies  of  confined  LO-phonon  modes  arc  determined 
for  a  GaAs/GaP  short-period  supcrluiiice  with  two  mono- 
layers  per  superlative  layer.  It  is  found  that  varying  inter¬ 
facial  force  constants  by  values  as  extreme  as  10%  results 
in  only  about  a  2%  change  in  the  frequencies'  of  the 
confined  LO-phonon  modes.  As  a  practical  matter, 
changes  in  the  frequencies  of  confined  LO-phonon  modes 
will  generally  be  considerably  less  than  2%  since  m  most 
supetlaiiiccs  and  quantum  wells  the  ratio  of  the  number  of 
bonds  at  interfaces  to  the  number  of  bonds  removed  trout 
interfaces  is  less  than  for  the  two-monolayer  per  layer  case 
examined  herein. 

2.  Linear-!  Ttait:  Model 
with  Variable  Inu-rfacial  I'Tirce  ( "rinslanls 

In  order  to  introduce  the  basic  properties  of  the 
linear-chain  model  used  in  this  analysis,  we  will  consider 
utili.illy  die  case  where  interface  ctfecis  are  absent. 
Figure  I  depicts  a  schematic  Inr  die  huear-ch.im  system 
used  10  study  the  role  of  foice-ciinwant  variations  at 
hctcro|unclioii  inter  faces.  The  z  axis  is  chosen  as  the 
growth  diiecliiin  ol  die  lattice  Hence,  for  the  longitudinal 
modes  of  interest,  the  atomic  displacements  aic  in  the  z 
directum  In  the  absence  of  strain,  the  force  constant  m 
each  layer  ts  estimated  based  on  the  frequencies  ot  die 
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Abstract  -  This  paper  describes  applications  which  combine  Monte  Carlo  methods 
with  other  techniques  to  model  semiconductor  devices.  These  procedures  can  produce 
more  effective  means  to  describe  the  behavior  of  device  structures  requiring  detailed 
physical  descriptions  and  increased  computational  efficiency.  Applications  presented 
in  this  paper  include:  (a)  simulation  of  transport  across  hetero-barriers  with  quasi- 
ballistic  effects;  (b)  simulation  of  metal-semiconductor-metal  photodetectors  in  which 
parasitic  circuit  effects  are  important;  and  (c)  generation  of  transport  parameters  for 
use  in  drift-diffusion  ( and  hydrodynamic)  models  while  negating  assumptions  about 
the  nature  of  the  particle  velocity  distribution  function. 

1.  Introduction 

The  Monte  Carlo  method  is  a  versatile  numerical  technique  for  studying 
semiconductor  materials  and  devices  (1-4).  It  provides  mechanisms  for  including 
detailed  physics  of  particle  transport  through  solid  state  materials,  especially  the 
clectron-phonpn  interaction  and  other  scattering  processes,  realistic  energy  band 
structure,  and  temporal-  and  spatial-dependent  parameters.  Steady-state  and  transient 
properties  arc  accessible  with  the  Monte  Carlo  method.  Realistic  boundary 
conditions  can  be  imposed  which  allow  simulation  of  a  variety  of  device  structures. 
In  principle,  when  coupled  with  a  self-consistent  numerical  solver  for  Poisson’s 
equation,  the  Monte  Carlo  method  permits  one  to  explore  die  properties  of  realistic 
two-  and  three-dimensional  semiconductor  devices  with  detail  and  precision. 
However,  as  more  physical  effects  arc  included,  as  spatial  dimensions  decrease  toward 
the  quantum  limit,  and  as  fluctuation  phenomena  become  more  important, 
limitations  to  the  Monte  Carlo  method  become  prohibitive  for  device  simulations. 
In  general,  this  limitation  first  manifests  itself  through  requirements  for  increased 
computational  resources  and  improved  numerical  algorithms. 

The  Monte  Carlo  method  generates  exact  numerical  solutions  to  the  Boltzmann 
transport  equation  by  stochastically  tracking  particle  movements  in  phase  space  (2). 
While  use  of  numerical  estimators  yield  parameters  of  physical  interest  (e.g.,  drift 
velocity),  the  fundamental  quantity  obtained  from  the  simulation  is  the  particle  phase 
space  distribution  function.  In  this  paper,  we  will  discuss  ways  in  which  distribution 
function  information  obtained  from  Monte  Carlo  simulations  can  be  used  to  augment 
the  study  of  device  behavior  using  other  simulation  methods. 

2.  Quasi-Ballistic  Transport  Across  Hetero-Barriers 

The  distribution  function  experiences  fundamental  changes  when  an  ensemble  of 
carriers  passes  dirough  an  abrupt  potential  barrier.  By  imposing  continuity  of  total 
particle  energy  and  the  parallel  component  of  wave  vector,  there  are  two  important 
consequences  on  the  velocity  distribution  function.  First,  the  normal  component  of 
the  velocity  vector  is  shifted  to  high  positive  values  by  an  jmount  equal  to  the 
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Abstract  -  This  paper  presents  a  new  hydrodynamic  transport  model  with  non¬ 
parabolic  conduction  bands  and  quantum  correction  terms.  For  the  first  lime 
solutions  for  the  full  quantum  balance  equations,  applied  to  an  ultrasmall  elec¬ 
tron  device,  are  presented. 

1.  Introduction 

Requirements  for  faster  electronics  have  produced  smaller  devices  influenced  by 
quantum  interference  effects  which  can  not  be  modeled  by  classical  theory.  These 
devices  operate  under  nonequilibrium  conditions  where  the  average  carrier  energy 
reaches  many  times  its  equilibrium  value.  To  address  these  problems  several  ap¬ 
proaches  [lj  have  been  considered.  This  paper  describes  a  preliminary  investigation 
into  solutions  for  electron  density,  average  velocity  and  average  energy  for  an  ultra¬ 
small  electronic  structure  under  the  conditions  of  nondassical  electron  transport, 
using  an  approach  previously  described  as  quantum  hydrodynamics  [2|. 

2.  The  Nottparnbolic  Model 

First,  we  present  a  brief  summary  of  a  unique  form  of  the  hydrodynamic  transport 
model  applicable  to  nonparabolic  conduction  bands.  A  complete  description  of  this 
model  will  be  formalized  elsewhere  [3).  Our  hydrodynamic  model  was  developed 
by  studying  moments  of  the  Boltzmann  transport  equation  (4|.  This  process  was 
achieved  using  the  moment  operators  4,|  —  1 .  'In  =  u(k)  =  "flcj'  anl*  '^3  “ 

I  rn(k)u(k)  ■  u(k)  /:i.(k).  Here,  the  Kane  dispersion  relation.  jj-S  n,-l  I  i 

a l'Je).  has  been  used  to  define  u(k)  and  >»i(k)  —  in’  \J  I  }-  These  moment 

operators  were  chosen  because  they  lead  to  a  form  which  can  be  manipulated 
more  easily  and  one  in  which  simplifying  approximations  can  be  seen  more  clearly 
Performing  the  moment  process  and  simplying  to  first  order  yields  immediately  the 
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Metal-semiconductor-metal  (MSM)  photoderectors  are  of  interest  for  integrated 
optoelectronic  circuits  (OEICs)  and  other  discrete  circuit  applications  because  this  very  high¬ 
speed  device  can  be  fabricated  using  conventional  GaAs  integrated  circuit  technology.  The 
typical  configuration  of  such  a  photodetector  consists  of  two  interdigitated  Schotiky  electrodes 
patterned  on  the  semiconductor  surface.  The  fastest  MSM  detector  has  an  intrinsic  bandwidth 
of  105  GHz.  which  compares  well  with  bandwidths  of  other  Schotiky  photodiodes.  The 
potential  to  achieve  even  larger  bandwidths  appears  to  be  quite  promising. 

We  have  performed  simulations  of  intrinsic  detectors  and  obtained  time  responses  for 
devices  with  interfinger  spacings  between  0.1-1  pm.  A  self-consistent  ensemble  Monte  Carlo 
simulation  method  was  employed  to  study  electron  and  hole  transport  under  low  illumination 
intensity.  The  Monte  Carlo  method  provides  mechanisms  for  including  detailed  physics  of 
carrier  transport,  especially  the  electron-phonon  interaction  and  other  scattering  processes, 
realistic  energy  band  structure,  and  temporal-  and  spatial-dependent  parameters.  Realistic 
boundary  conditions  are  imposed  and  the  utilization  of  a  self-consistent  numerical  solver  for 
Poisson's  equation  permits  us  to  simulate  realistic  one-  and  two-dimensional  structures  with 
detail  and  precision.  In  addition,  the  time-dependent  Schrodinger  equation  is  solved  for  the 
detector  with  0.1  pm  distance  between  fingers  under  an  assumption  of  collisionless  carrier 
transport,  and  a  comparison  between  the  two  different  methods  is  made.  The  devices  are 
dominated  by  stationary  hot  electron  and  hole  transport. 

A  parasitic  circuit  is  introduced  to  more  closely  model  realistic  devices.  Tlte  intrinsic 
responses  are  used  in  a  circuit  simulator  to  obtain  the  time  responses  of  pliotodetectors  due  to 
excitation  by  a  very  short  optical  pulse.  The  spectral  power  density  is  calculated  using  a  fast 
Fourier  transform  (FFT)  to  predict  the  bandwidth  of  intrinsic  and  extrinsic  photodctcctors. 
This  simulation  technique  is  shown  to  yield  good  agreement  with  published  experimental  data 
for  a  GaAs  photodetector. 

For  GaAs  MSM  photodctcctors,  we  predict  that  bandwidths  in  excess  of  150  GHz  can  be 
achieved  for  devices  with  0.1 -0.2  pm  interfingcr  spacing.  The  parasitic  effects  of  the  contact 
pattern  arc  critical  in  the  performance  of  such  short  devices,  and  further  optimization  of  the 
interfingcr  coplanar  transmission  line  could  result  in  even  higher  bandwidths.  The  role  of  the 
transient  response  of  holes  on  the  intrinsic  and  extrinsic  photodetector  response  and  bandwidth 
will  be  discussed.  Results  arc  also  presented  for  InAlAs/InGaAs/inP  MSM  photodctcctors. 
These  devices  operate  at  longer  wavelengths  than  GaAs  MSM  photodctcctors.  and  are  thought 
to  provide  potential  for  high-speed  near-infrared  transceivers.  Our  results  predict  that  hole 
transport  and  transient  response  limit  the  bandwidth  achievable  with  this  material  system  to 
about  25  GHz  for  MSM  photodctcctors  with  0.2  pm  interfingcr  spacing  and  current  contact 
design  geometries. 
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The  polarization  eigenvectors  and  dispersion  relations  for  confined  longitudinal-optical  (LO) 
and  surface-optical  (SO)  phonon  modes  in  a  rectangular  quantum  wire  are  derived  within  the 
continuum  approximation.  The  SO  and  confined  LO  modes  are  compared  for  a  quantum  wire 
with  a  square  cross  section  in  the  limit  of  a  vanishing  phonon  wave  vector  along  the  quantum-wire 
axis.  For  this  case,  it  is  demonstrated  that  the  SO  mode  becomes  dominant  over  the  confined  LO 
modes  in  interaction  with  electrons  and  must  be  taken  into  account  in  calculating  the 
electron-optical-phonon  scattering  rate. 


Epitaxial  techniques  for  the  growth  of  compound-  cently  have  confined  LO-phonon  modes  been  studied  in 

semiconductor  structures  have  advanced  to  the  level  more  complex  structures  such  as  quantum  wires4  and 

where  it  is  possible  to  fabricate  wire-like  regions  of  low-  strained-layer  short-period  superlatticcs.10 11  In  this 

band-gap  semiconductor  material  surrounded  completely  study,  the  polarization  eigenvectors  for  the  confined  LO- 

by  regions  of  higher-band-gap  semiconductor  material.  In  phonon  and  the  SO-phonon  modes  of  a  quantum  wire  sur- 

particular,  such  wire-like  structures  have  been  realized  rounded  by  a  medium  with  unity  dielectric  constant  are 

with  rectangular  cross  sections  having  small  dimensions  derived  in  the  continuum  approximation.4  5  7-9  For  pho- 

relative  U.  the  thermal  de  Broglie  wavelength.1  In  these  non  wavelengths  long  compared  with  the  lattice  constant, 

structures  the  electron-optical -phonon  scattering  rate  is  continuum  models  are  expected  to  be  valid5 ' 9  and  the  rc- 

affectcd  not  only  by  changes  in  the  c'cclron  wave  function  suits  of  such  models  agree  well  with  those  of  microscopic 

due  to  the  confining  rectangular  potential  but  also  by  models.'2  The  quantum  wire  is  taken  to  have  a  dielectric 

changes  in  the  longitudinal-optica!  'LO)  phonon  modes  constant  and  to  be  bounded  by  ±  Lr! 2  and  ±  Ljl  in 

caused  by  phonon  confinement.  Fas-,  ei  air  have  recently  the  y  and  z  directions,  respectively.  The  region  surround- 

presented  striking  experimental  evidence  of  phonon  ing  the  quantum  wire  is  taken  to  have  a  dielectric  constant 

confinement.  Size  effects  on  the  total  scattering  rates  for  e,  (*“l).  The  confined  LO-phonon  modes  for  such  a 

polar-optical -phonon  scattering  of  onc-dimensiona!  (ID)  quantum  wire  have  been  given  previously.4 
electron  gases  .n  quantum  wires  and  two-dimensional  Since  this  system  is  translational!}  invariant  in  the  x 
(2D)  clcaron  gases  In  quantum  wells  have  been  evaluated  direction,  the  potential  describing  the  optical  phonon 
previously  by  Lcburton. J  Recently,  Leburton's  treatment  modes  may  be  taken  as 
of  ID  clcctron-LO-phonon  scattering  was  extended  by  re¬ 
placing  bulk  LO-phonon  modes  with  the  confined  LO-  <J>(r)  — <J>(>\:)e,i,x .  (l) 

modes  of  a  quantum  wire.4  The  scattering  rate  calcula¬ 
tions  of  Ref.  4  included  confined  LO-phonon  modes  of  a  where  is  the  phonon  wave  vector  in  the  x  direction.  In 

rectangular  wire  and  neglected  electron  scattering  due  to  thc  abscncc  of  frec  chargC(  {hc  dlv.crgcnce  0f  ihc  dis. 

the  surface-optical  (SO)  phonon  modes  at  the  quantum-  placement  vector  must  vanish  and  it  follows  that  the  po- 

wire  boundaries.  In  this  Rapid  Communication,  the  lcnliaI  «x>(y,r)  of  thc  phonon  modes  must  satisfy 

confined  LO-  and  SO-phonon  modes  are  derived  for  a  rec¬ 
tangular  quantum  wire.  In  thc  limit  of  small  (phonon) 
wave  vectors  along  the  quantum-wire  axis,  it  is  shown  that  b‘  ^  d: 

thc  contribution  of  clcctron-SO-phonon  scattering  to  the  Qy'  q-' 

ground-state  clcctron-LO-phonon  scattering  rate  be¬ 
comes  significant  and  cannot  be  neglected  compared  to 
the  contribution  by  the  confined  LO-phonon  modes.  where  z»  for  -  Ljl  <y  <+/.,.  "2  and  - l.-/2< : 

Since  ihc  early  work  of  Fuchs  and  Klicwcr.56  thc  <  •*•/„./ 2.  outside  of  this  region  whiih  defines  thc  qu„.. 
effects  of  confinement  on  LO-phonon  .nodes  have  been  turn  wire,  e  **«,.  Thc  derivation  of  Lq  (2)  is  straight!  a 
studied  theoretically  by  a  number  of  authors  For  the  ease  wjrd  upon  taking  E(r)  "  -  V«I»(r).  Dtr)  ^iEU)  -  E'r 
of  phonon  confinement  in  a  quantum  well.’ Onlj  re-  f  4.7P(r),  ar.d  «I»fr)  as  given  by  Lq.  » !).  ..  .unil.u  u-.a 
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The  roles  of  DX  centers  and  surface  states  associated  with  the  //-AlGaAs  layer  of  fi-doped 
AlGaAs/GaAs  high-elect ron-mobility  transistors  have  bee-  investigated  by  employing 
a  two-dimensional,  self-consistent  ensemble  Monte  r—to  simulation.  It  is  found  that  both  DX 
centers  ail'5  surface  states  degrade  device  performance,  particularly  as  gate-to-sourcc 
voltage  increases.  This  degradation  is  manifested  largely  through  reduced  channel  electron 
concentration,  increased  intcrvalley  scattering,  and  enhanced  real-space  transfer.  Of 
the  two  defect  states,  DX  centers  ha  more  detrimental  influence  on  drain  current  and 
transconductance. 


I.  INTRODUCTION 

Recent  studies1"1  of  the  AlGaAs/GaAs  high-electron- 
mobility  transistor  (HEMT)  have  indicated  that  deep  do¬ 
nor  levels  (DX  centers)  and  surface  states  associated  with 
the  //-AlGaAs  layer  can  degrade  device  performance  sig¬ 
nificantly.  Experimental  measure;'  us  show  that  there  is  a 
sharp  increase  in  the  concentration  of  trapped  carriers  by 
DX  centers  (and  a  subsequent  collapse  in  device  transcon¬ 
ductance)  in  (he  AlAGa|  _  aAs  layer  with  Al  mole  fraction 
*>0.2  Similar  degradation  has  been  observed  in  connec¬ 
tion  whli  surface  states  as  well.  The  nature  of  DX  centers 
:  i  AlGaAs  has  been  a  subject  of  controversy,  investigated 
by  many  authors  experimentally.5'7  Mizuta  et  al ,5  describe 
the  DX  center  as  a  simple  substitutional  donor,  while  the 
work  by  Dhar  et  aib  reveals  that  there  are  two  closely 
spaced  deep  donor  levels  in  addition  to  a  shallow  donor 
level  Study  also  shows  tins'  the  ratio  of  shallow  donor 
density  to  deep  donor  density  exhibits  substantial  varia¬ 
tions  from  sample  to  sample,  depending  upon  the  material 
composition  and  (possibly)  some  other  g.owth  parame¬ 
ters.1 

To  date,  theoretical  hives' igations4,8'7  of  the  effects  of 
DX  centers  and  surface  stales  on  device  performance  have 
been  scarce  and  largely  inconclusive.  In  particular,  most  of 
the  efforts  have  been  devoted  to  the  phenomenological  de 
scription  (i  e  .  parameter  fitting)  of  device  degiadation  in 
the  presence  of  eithei  DX  centers  or  surface  states,  rather 
than  to  physical  models  which  can  elucidate  the  underlying 
mechanisms  The  combined  trapping  effects  by  both  DX 
centers  and  surface  states  have  n  l  been  identified  and,  in 
addition,  are  not  well  understood. 

Ii  is  the  purpose  of  this  paper  to  present  lcsalls  of  our 
study  on  the  influence  of  DX  centers  and  su  fare  states 
employing  a  realistic,  physical  model.  The  5  doped 
AlGaAs/GaAs  MKMT  has  bee.i  chosen  as  a  specific  struc¬ 
ture  to  be  studied  due  to  its  superior  performance  and 
recent  interest  compared  to  lire  uniformly  doped 
counterpart  1,1  11  'fire  collective  roles  of  DX  centers  anil 
surface  states  arc  examined  for  0.3  /on  gale-length  fi-dop.d 
AIGaAs/( J..As  HEM  I  s  with  gerreial  napping  situations 


and  different  bias  conditions  by  employing  a  two- 
dimensional,  self-consistent  ensemble  Monte  Carlo  simula¬ 
tion.  The  results  clearly  indicate  the  importance  of  DX 
centers  and  surface  states  on  device  modeling  and  provide 
an  explanation  for  detailed  degiadation  mechanisms. 

II.  MODEL 

The  Monte  Carlo  model  incorporates  three-valley  (T- 
L-X)  nonparabolic  analytical  band  structures  for  both 
AlGaAs  and  GuAs.  The  Irajcetories  of  about  I.SxlO4 
sample  electrons  are  traced  both  in  real  space  and  in  mo¬ 
mentum  space  under  the  influence  of  a  two-dimensional 
electric  field  generated  by  a  self-consistent  Poisson  equa¬ 
tion  solver.  These  electrons  arc  subject  to  relevant  scatter¬ 
ing  mechanisms  such  as  polar  optical,  intervalley,  lomzed- 
impurity,  and  electron-electron  scatterings.  The 
contribution  by  the  collective  excitations  of  the  electron 
gas  (i.e.,  plusmou)  is  considered  through  frequent  updat¬ 
ing  of  (lie  electric-field  profile  (every  2x10  15  s).  Real- 
space  transfer  of  electrons  between  GaAs  and  AlGaAs  is 
included.  Since  the  two-dimensional  quantum  effects  are 
relatively  unimportant  for  the  electrons  in  the  high-field 
region,14  we  treat  the  AlGaAs/GaAs  interface  potential 
classically  and  ignore  size  quantization.  The  effect  of  cai- 
lier  degeneracy  is  taken  into  account  for  elections  in  the  f 
valley  by  employ  mb  Fermi  statistics  with  the  calculated 
local  quasi-Fcimi  level  and  F -valley  average  election  en- 
cigy.l,!5  Chaigc* neutrality  in  c.  small  region  adjacent  to 
the  source  and  d’ain  contacts  is  maintained  lluoughout  the 
simulation  and  serves  as  the  cuteiion  foi  ejection  injection. 
Figure  1  shows  a  schematic  drawing  of  typical  ().5-/im 
t V  doped  AlGaAs/GaAs  HEM  I  structure  studied  in  this 
paper,  'flic  device  consists  of  a  -400- A  unintentionally 
doped  (5  a  1014  cm  ')  top  AlGaAs  layei,  a  UK)- A  xlcp- 
doped  (5  x  1018  cm  J)  //  *  -AlGaAs  layei,  and  a  50  A 
uniiitcnlionally  doped  (5  a  I0n  cm  ')  AlGaAs  spacer 
y;i  followed  by  a  2()ot)-A  unintentionally  doped  (5 
a  Id14  cm  ')  GaAs  layei.  The  contact  doping  is  assumed 
to  be  I  a  l(in  cm  '.  An  Al  .nole  had  ion  of  0.3  is  used 
loi  each  AlGaAs  layei.  Hie  gntc-clcctiodc  Sehollky  bai- 
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ABSTRACT 

The  Hamiltonian  describing  the  interaction  of  both  confined  longitudinal-optical  and  surface-optical 
phonons  with  charge  carriers  is  derived  from  the  macroscopic  dielectric  continuum  model  for  the  case  of 
a  rectangular  quantum  wire  where  phonon  confinement  occurs  in  two  of  the  three  spatial  dimensions. 
The  full  interaction  Hamiltonian  is  used  to  calculate  the  total  scattering  rate  for  electron-optical-phonon 
scattering  of  electrons  traversing  a  GaAs  square  quantum  wire.  The  results  demonstrate  that  the  interac¬ 
tion  by  the  surface-optical  phonon  modes  is  very  strong  and  may  dominate  over  other  scattering 
processes,  especially  with  dimensions  of  about  100  A  or  less.  A  considerable  decrease  in  the  total 
scattering  rate  for  optical  phonons  as  a  result  of  simple  reduction  in  dimensionality  is  not  observed  in  this 
study. 


1.  INTRODUCTION 


Epitaxial  techniques  for  the  growth  of  compound-semiconductor  structures  have  advanced  to  the 
level  where  it  is  possible  to  fabricate  wire-like  regions  of  narrow-band-gap  semiconductor  material  sur¬ 
rounded  completely  by  regions  of  large-band-gap  semiconductor  material,  in  particular,  such  wire-like 
structures  have  been  fabricated  with  rectangular  cross  sections  having  dimensions  small  relative  to  the 
electron  thermal  de  Broglie  wavelength;  furthermore,  optical  anisotropy  measurements  on  these  wire-like 
structures  provide  evidence  of  two-dimensional  quantum  confinement  of  charge  carriers.1  In  such  struc¬ 
tures,  the  electron-optical-phonon  scattering  rate  is  affected  by  changes  in  the  wave  functions  of  the 
charge  carriers  due  to  the  confining  rectangular  potential  but  also  by  changes  in  the  carrier-optical-phonon 
interaction  Hamiltonian  caused  by  phonon  confinement  and  localization  of  charge  carriers  due  to  the 
confining  rectangular  potential.  Indeed,  recent  measurements  on  cylindrical  quantum  wires2  have  pro¬ 
vided  striking  evidence  of  surface-optical  (SO)  modes  and  Fasol  et  al.3  have  experimentally  verified  pho¬ 
non  confinement.  The  effects  of  carrier  confinement  on  the  total  scattering  rates  for  poiar-optical-phonon 
scattering  in  quantum  wells  and  quantum  wires  have  been  evaluated  previously  by  Leburton4  with  bulk 
longitudinal-optical  (LO)  phonon  modes. 

In  this  paper,  one-dimensional  (ID)  clectron-opiical-phonon  scattering  rates  in  rectangular  quantum 
wires  are  determined  by  including  the  effects  of  carrier  confinement  as  well  as  by  replacing  the  bulk  LO- 
phonon  modes  by  the  confined  LO-phonon  modes  and  s'1  ^ace-optical  (SO)  phonon  modes  for  a  rectangu¬ 
lar  quantum  wire.  While  optical-phonon  confineme’  •  eti  cts  are  clearly  evident  in  quantum-well  struc¬ 
tures,3  the  results  of  this  paper  demonstrate  (hat  can.,--'-  ">-phonon  scattering  may  dominate  over  other 
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Abstract— Transport  properties  and  device  performance  of 
delta-doped  and  uniformly  doped  AIGaAs/GaAs  high  electron 
mobility  transistors  (HEMT’s)  with  identical  threshold  voltages 
and  gate  capacitors  are  investigated  with  two-dimensional,  self- 
consistent  ensemble  Monte  Carlo  simulations.  The  model  in¬ 
cludes  the  effects  of  real-space  transfer  and  carrier  degeneracy, 
as  well  as  the  influence  of  DX  centers  and  surface  states.  A  one- 
to-one  comparison  of  simulation  results  for  the  two  devices 
demonstrates  superior  performance  for  the  delta-doped  HEMT 
and  provides  physical  basis  for  the  observed  improvements.  In 
particular,  the  delta-doped  HEMT  maintains  its  superior  de¬ 
vice  performance  as  gate  bias  is  increased.  The  simulations 
provide  reasons  for  these  improvements.  The  predicted  advan¬ 
tages  of  the  delta-doped  HEMT  include  high  channel  electron 
density  and  high  channel  drift  velocity,  which  stem  primarily 
from  improved  channel  electron  confinement,  reduced  parallel 
conduction  at  high  gate-to-source  voltage,  and  reduced  peak 
channel  electric  field.  These  advantages  of  delta-doped  HEMT 
are  manifested  in  higher  transconductance  and  improved  drain 
current  drive  capability  at  a  wider  gate-to-source  voltage,  which 
are  desirable  properties  for  ultra-fast  digital  and  microwave 
device  applications. 

I.  Introduction 

HE  HIGH  electron  mobility  transistor  (HEMT)  is  a 
promising  candidate  for  ultra-fast  digital  and  micro- 
wave  device  applications.  Conventional  uniformly  doped 
AIGaAs/GaAs  HEMT’s  are  superior  to  comparable 
MESFET  structures  in  overall  device  performance.  How¬ 
ever,  further  device  improvement  is  limited  by  the  occur¬ 
rence  of  persistent  photocondutivity  [1],  threshold-volt¬ 
age  shift  [2],  and  collapse  of  l-V  characteristics  [3]  due 
to  effects  of  deep  donors  in  the  AlGaAs  layer.  The  utili¬ 
zation  of  low  A1  mole  fraction  can  reduce  these  effects; 
however,  this  also  reduces  the  conduction-band-edge  dis¬ 
continuity,  resulting  in  decreased  channel  electron  den¬ 
sity.  This  inherent  limitation  can  be  circumvented  either 
by  device  structure  modifications  (such  as  delta-doping 
[4]-[8])  or  by  using  pseudomorphic  heterostructures  (9), 
[10].  Recently,  much  attention  has  been  paid  to  delta- 
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doped  (also  referred  to  as  pulse-doped  [4],  [5]  or  planar- 
doped  [7])  HEMT  structures  due  to  the  maturity  of  epi¬ 
taxial  technology  and  excellent  lattice  match  for  the 
AIGaAs/GaAs  sj  ,tem.  By  introducing  an  unintentionally 
doped  i-laycr  between  Schottky  gate  and  n* -AlGaAs 
layer,  the  delta-doped  HEMT’s  provide  high  channel 
electron  density,  reduced  trapping  effect,  and  improved 
threshold  voltage  control  as  well  as  high  breakdown  char¬ 
acteristics.  Also,  one  expects  improved  current  drive  ca¬ 
pability  and  high  transconductance  due  to  high  electron 
density  and  small  gate-to-channel  spacing  for  such  de¬ 
vices. 

Although  the  advantages  of  the  delta-doped  HEMT  over 
the  uniformly  doped  HEMT  have  Men  shown  experimen¬ 
tally,  detailed  theoretical  invest  ..ions  are  essential  to 
fully  understand  underlying  physics  and  exploit  advan¬ 
tages  of  delta-doping  on  electron  transport  and  device  per¬ 
formance  in  the  submicrometer  regime.  While  systematic 
analysis  of  conventional  uniformly  doped  AIGaAs/GaAs 
HEMT’s  and  pseudomorphic  HEMT’s  have  been  re¬ 
ported  by  many  authors  [  1 1]— [  14],  such  studies  on  delta- 
doped  structures  have  been  rare  and  mainly  first-order  [7], 
[8],  It  is  well  known  that  the  operation  of  submicrometer 
heterostructure  devices  involves  highly  nonlinear  trans¬ 
port,  such  as  hot-electron  effects  and  real-space  transfer, 
and  requires  a  powerful  device  model  to  accurately  de¬ 
scribe  carrier  transport  behavior  and  device  performance. 
One  of  the  most  sophisticated  and  complete  tools  is  the 
ensemble  Monte  Carlo  simulation  coupled  with  a  self- 
consistent  Poisson  solver.  In  this  study,  we  examine  the 
transport  properties  and  device  performance  of  compara¬ 
ble  0.5-fim  gate-length  delta-doped  HEMT  and  uniformly 
doped  HEMT  (hereafter,  D-HEMT  and  U-HEMT,  re¬ 
spectively)  using  a  self-consistent  two-dimensional  en¬ 
semble  Monte  Carlo  simulation.  Our  simulation  results 
demonstrate  the  superior  performances  of  D-HEMT’s  to 
those  of  comparable  U-HEMT’s  in  terms  of  device  struc¬ 
ture,  doping,  and  bias  conditions  and  provide  a  clear 
physical  explanation  for  the  improvements. 

II.  Model 

The  AIGaAs/GaAs  HEMT  structure  used  in  the  sim¬ 
ulations  is  illustrated  in  Fig.  I.  For  th  )-HEMT.  the 
device  consists  of  a  400-A  unintentionally  doped  (5  x 
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The  Feynman  influence  functional  for  the  polar  coupling  of  carriers  to  confined  phonon  modes  in  het- 
erostructurcs  is  considered.  It  is  found  that  the  partial  contributions  to  the  total  influence  functional 
from  the  individual  branches  of  the  optical-phonon  energy  spectrum  are  related  and  constrained  by  con¬ 
servation  relations.  The  derivation  of  these  conservation  relations  requires  no  specific  functional  form 
for  the  phonon  modes;  rather  it  employs  only  the  inherent  orthogonality  and  mathematical  completeness 
of  the  classical  vibrational  modes  over  the  crystal-lattice  degrees  of  freedom.  In  the  bulk-crystal  limit 
these  conservation  relations  lead  to  the  familiar  results  of  Feynman;  for  arbitrary  heterostructures,  these 
conservation  relations  provide  a  basis  for  an  estimation  of  the  influence  functional  from  limited 
knowledge  of  the  spatial  confinement  of  the  individual  phonon  branches. 


I.  INTRODUCTION 

The  Feynman  path-integral  (FPI)  formulation  of  quan¬ 
tum  mechanics  allows  formal  inclusion  of  carrier-phonon 
coupling  to  all  orders  in  the  coupling  potential  via  the 
influence  functional.1,2  For  this  reason,  the  FPI  formal¬ 
ism  long  has  been  used  to  calculate  carrier  self-energies 
and  effective  masses  in  bulk  semiconductors.3"6  More 
recently,  the  FPI  formalism  has  been  extended  to  the 
study  of  confined  carriers  in  ultrasmall  heterostructure 
geometries,7  including  recent  advances  in  numerical 
methods.3"10  However,  in  heterostructures,  beyond  the 
confinement  of  carriers,  the  confinement  of  phonons 
should  be  considered  as  well.  Here,  this  latter 
confinement  effect  on  the  influence  functionals  for  polar 
coupling  of  carriers  to  longitudinal-optical  (LO)  phonons 
is  considered,  without  assuming  any  specific  functional 
form  for  the  phonon  modes.  In  Sec.  II,  the  influence 
functional  is  defined  in  terms  of  the  carrier-phonon  in¬ 
teraction  Hamiltonian.  In  Sec.  Ill,  it  is  shown  that  the 
partial  contributions  to  the  total  influence  functional 
from  the  individual  branches  of  the  LO-phonon  spectrum 
are  related  and  constrained  by  conservation  relations.  In 
Sec.  IV,  it  is  shown  that  these  conservation  relations  lead 

_ ! 


to  tne  familiar  results  of  Feynman3,4  in  the  bulk-crystal 
limit,  and  that  for  heterostructures  they  provide  a  basis 
for  estimation  of  the  influence  functional  from  limited 
knowledge  of  the  spatial  confinement  of  the  individual 
phonon  branches. 

II.  PATH-INTEGRAL  FORMALISM 


Consider  a  coupled  carrier-phonon  system  with 
quantum-mechanical  Hamiltonian  of  the  form 

ft1  ■> 

H= - - — V2+  F(r) 

2m*  r 


+  2 


2  di 


(1) 


where  r  is  the  carrier  coordinate,  tfm  are  the  phonon 
coordinates,  ftoom  are  the  uncoupled  phonon  energies,  and 
<im(r)  are  real  functions  characterizing  the  carrier  posi¬ 
tion  dependence  in  the  carrier-phonon  interaction.  For 
this  system,  the  equilibrium  density  matrix  for  the  carrier 
averaged  over  the  phonon  coordinates  and  subject  to  a 
Maxwell-Boitzmann  distribution  in  energies  at  tempera¬ 
ture  7*  is  given  by 


p(r,r'\T)  =  f'cxp 

J  r’ 


1  .7’ 

(4-w  *|r(/),2  +  V{r(t))]dt 

ft  J  o 


7{t(i)}7)tU) 


•  2) 


in  the  FPI  formalism.2  Here,  “  f' .  •  ■  ■,'/jrU)’’  signifies  an  integral  over  all  paths  from  r'  to  r  that  are  continuous  (unc¬ 
tions  of  time,  and  ;7[r(/)]  is  the  influence  functional.  7[r(/))  accounts  for  coupling  between  the  carriers  and  phonons  to 
all  orders  m  the  coupling  potential  and  is  given  by4 
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1  Introduction 


The  Monte  Carlo  method  is  a  powerful  technique  for  investigating  electron 
transport  physics  in  semiconductors  and  semiconductor  devices[l,  2,  3].  Monte 
Carlo  is  a  general  statistical  numerical  method  for  solving  mathematical  and 
physical  problems  which  was  firmly  established  long  before  being  introduced 
to  semiconductor  physics[4).  In  the  case  of  electron  transport,  the  method 
provides  the  most  efficiently  implemented  and  physically  adaptive  technique 
currently  available  to  determine  exact  numerical  solutions  of  the  electron  motion 
in  semiconducting  materials  and  devices[5,  and  references  therein]. 

The  Monte  Carlo  method  provides  a  relatively  simple  and  accurate  indirect  ap¬ 
proach  to  determine  particle  distribution  functions  and  transport  parameters 
resulting  from  the  Boltsmann  tranport  equation  subject  to  complex  scatter¬ 
ing  mechanisms[6,  7).  In  the  method,  a  type  of  mathematical  experiment  is 
performed  where  the  trajectory  of  an  electron  (or  an  ensemble  of  electrons)  is 
simulated  subject  to  applied  fields  and  given  scattering  mechanisms.  The  du¬ 
ration  of  the  carrier  free  flights  and  the  particular  collision  events  are  selected 
stochastically  (randomly)  according  to  known  scattering  probabilities  which  de¬ 
scribe  the  microscopic  processes.  When  an  electron  is  tracked  over  a  tufficienlly 
long  history  the  dynamics  can  then  be  used  to  estimate  lime-independent  prop¬ 
erties  of  the  entire  electron  gas  (For  time  and/or  space  dependent  problems, 
the  sample  ensemble  must  be  sufficiently  large  to  accurately  represent  the  en¬ 
tire  electron  gas).  This  procedure  permits  the  investigation  and  interpretation 
of  physical  phenomena  unattainable  solely  by  experimental  methods.  Thus,  the 
Monte  Carlo  method  is  extremely  valuable  to  solid  stale  physics  from  both  a 
theoretical  and  experimental  point  of  view. 
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Dielectric  continuum  models  of  optical-phonon  modes  predict  an  enhancement  in  the 
strength  of  the  surface-optical  (SO)  modes  in  double-barrier  hetcrostructurcs  as  the 
hcterojunction-to-hetcrojunction  separation  is  reduced.  There  is  currently  no 
consensus  on  the  nature  of  the  elect ron-SO-phonon  coupling  interaction.  In  this  work,  the 
ratio  of  electron  scattering  by  the  SO-phonon  modes  to  that  by  the  confined 
longitudinal-optical  (LO)  phonon  modes  is  calculated  for  a  GaAs/AlAs  short-period 
superlatticc  based  on  the  assumption  that  the  electron-SO-phonon  interaction  may  be  described 
by  a  scalar  potential.  The  scaling  of  the  ratio  of  electron-SO-phonon  scattering  to  elect  ron- 
LO-phonon  scattering  as  a  function  of  the  superlatticc  period  provides  a  sensitive 
test  of  the  appropriateness  of  the  scalar-potential  model. 


Electron  interactions  with  longitudinal-optical  (LO) 
phonon  modes  in  hetcrostructurcs  arc  strongly  affected  by 
the  changes  in  the  Frohlich  Hamiltonian  caused  by  phonon 
confinement  and  localization,  as  well  as  by  the  changes  in 
the  electronic  wave  function  due  to  the  confining  potential. 
The  presence  of  hctcrointcrfaccs  gives  rise  to  the  confine¬ 
ment  of  optical  phonons  in  each  layer  (i.c.,  confined  mode) 
and  the  localization  in  the  vicinity  of  interfaces  [i.c.,  inter¬ 
face  mode  or  surface-optical  (SO)  mode].  There  have  been 
suggestions  that  interactions  by  the  interface  modes  can  be 
significant  and  that  the  scattering  rate  due  to  the  confined 
modes  can  be  considerably  reduced  in  some  structures 
compared  to  the  bulk  LO-phonon  scattering  rate.1'3  There¬ 
fore,  an  appropriate  treatment  of  the  optical-phonon 
modes  in  quantized  systems  is  essential  for  the  understand¬ 
ing  of  electron  transport  in  hetcrostructurcs.  Recently, 
both  macroscopic  and  microscopic  approaches  to  elect  ron- 
optical-phonon  interactions  in  hetcrostructurcs1'"  have 
been  applied  in  theoretical  treatments.  Enhanced  electron- 
SO-phonon  (i.e.,  clcctron-inlerfacc-phonon)  scattering  in 
polar  semiconductors  with  confining  dimensions  less  than 
about  50  A  has  been  indicated  recently12,13  on  the  basis  of 
a  scalar-potential  modes.  The  appropriateness  of  using  a 
scalar  potential  to  model  the  clcctron-SO-phonon  interac¬ 
tion  has  been  questioned."  In  this  letter,  we  model  the 
clcctron-SO-phonon  interaction  with  a  scalar  potential  to 
calculate  the  ratio  of  clcctron-SO-phonon  scattering  to 
clcctron-confined-LO-phonon  scattering  in  a  short-period 
GaAs/AlAs  superlatticc.  The  results  provide  a  criterion 
which,  in  conjunction  with  experiments,  can  be  used  to 
examine  the  validity  of  the  scalar-potential  model. 

A  GaAs/AlAs  short-period  superlatlice  is  considered 


in  this  study.  Each  GaAs  layer  has  width  a  and  each  AlAs 
has  width  b.  We  will  consider  the  case  where  a  uniform 
electron  current  J  is  incident  parallel  to  the  superlatticc 
hctcrojunctions.  The  ratio  of  electron-SO-phonon  scatter¬ 
ing  to  clectron-LO-phonon  scattering  in  the  GaAs/AlAs 
superlatticc  is  calculated  through  the  Fermi  golden  rule  by 
treating  the  clectron-LO-phonon  and  clcctron-SO-phonon 
interaction  Hamiltonians4,6  as  perturbation  Hamiltonians 
and  by  taking  (he  electronic  wave  function  in  the  z  direc¬ 
tion  as  that  given  by  Cho  and  Prucnal  for  the  maximum 
edge  of  the  first  subband;15  that  is, 


where  k(!  is  the  electron  wave  vector  parallel  to  the  super¬ 
lattice  hctcrojunctions,  L 2  is  the  hetcrojunction  area  (un¬ 
der  the  assumption  that  L^a  -)•  b  and  L—  ).  r  defines  the 
parallel  components  of  the  position  vectors  in  the  supcrlat- 
ticcs,  and  the  other  quantities  arc  as  discussed  above.  The 
origin  of  the  z  axis  (i.c.,  growth  direction)  is  chosen  at  the 
middle  of  a  GaAs  layer.  Taking  the  electron-optical-pho- 
non  interaction  Hamiltonians  to  be  as  in  Ref.  6.  the  matrix 
elements  for  the  emission  of  the  lowest  order  confined  LO 
phonons  and  the  symmetric  phonons  arc  given,  respec¬ 
tively,  by 

Mas.  >  =  <!<'.  <V,  +  I  |//c(.y .  ,!k..V„)-  (2) 

where  /Vq  is  the  phonon  occupation  number  and  k’  is  the 
electron  wave  vector  after  phonon  emission.  Performing 
the  indicated  overlap  integrals,  Mt  becomes 
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